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Thermische Beschreibung der Hadronenproduktion in e+e 
KoUisionen 

Die vorliegende Arbeit beschaftigt sich mit der Hadronenproduktion in e^e" KoUisionen 
bei verschiedenen Schwerpunktsenergien im Rahmen des Statistischen Models eines Ha- 
dronen Resonanz Gases. Das Model, welches fiir diese Analyse verwendet wird, basiert 
auf dem kanonischen Ensemble mit der exakten Erhaltung von fiinf Quantenzahlen. Die 
entsprechende kanonische Zustandsgleichung wird in Quantenstatistik hergeleitet und der 
kanonische Faktor wird, als eines der grundlegenden Eigenschaften der kanonischen Be- 
schreibung, detailiert untersucht. Die Parameter des Models werden durch einen x^-Fit 
an die neuesten gemessenen Multiplizitaten bei ^/s=\0, 29-35, 91 und 130-200 GeV ge- 
wonnen. Wie in SchwerionenkoUisionen werden die Messwerte mit der reinen thermischen 
Produktion verglichen. In einem zweiten Szenario werden harte KoUisionen beriicksichtigt, 
die eine Verstarkung der Produktion von schweren Quarks zur Folge haben. Die Ergebnis- 
se zeigen, dass innerhalb der hohen Genauigkeit der Experimente von wenigen Prozent, 
keines der Datensatze zufriedenstellend mit dem gewahlten Ansatz erklart werden kann. 
Deshalb ist es fraglich, ob die Teilchen Produktion in e+e^ KoUisionen thermischen Ur- 
sprung hat. AUerdings wurde auch beobachtet, dass die Messwerte fur charm und bottom 
Teilchen, die in harten KoUisionen erzeugt werden durch grosse kanonische Verstarkung 
solcher Teilchen in schweren Jets nahezu unabhangig von den Modell Parametem sehr gut 
erklart werden konnen. 

Thermal description of hadron production in e+e~ colli- 
sions 

This thesis gives a comprehensive analysis of hadron production in e+e~ collisions at dif- 
ferent center of mass energies in the framework of the Statistical Model of the hadron 
resonance gas. The model used for the analysis is formulated in the canonical ensemble 
with exact conservation of five quantum numbers. The corresponding canonical partition 
function in quantum statistics is derived and the canonical factor as the striking feature of 
the canonical framework will be investigated in detail. The parameters of the underlying 
model were determined using a fit to the average multiplicities of the latest measurements 
at y^=10, 29-35, 91 and 130-200 GeV. The measurements are compared to the pure ther- 
mal production like in heavy ion collisions as well as to the production in a second scenario 
accounting for hard collisions and the resulting enhancement of heavy quark. The results 
demonstrate that, within the accuracy of the experiments, none of the data sets is satisfacto- 
rily described with this approach, calling into question the notion that particle production 
in e"'"e~ collisions is thermal in origin. Nevertheless it is observed that the charm and 
bottom particle yields produced in hard collisions can be explained by large canonical en- 
hancements of such particles in heavy jets, almost independent of the model parameters. 
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Chapter 1 
Introduction 



Quantum Chromodynamics (QCD) is a non-Abelian gauge theory describing the strong 
interaction which takes place between color charged particles called quarks and gluons. 
The strong coupling constant is by far not a real constant, but depends on energy. Unlike 
the electromagnetic coupling constant in QED the strong coupling constant is increasing 
with decreasing energy. Therefore the strong interaction has some special features like 
confinement and asymptotic freedom. Asymptotic freedom means that as distance scales 
decrease and energy scales increase the coupling constant for the strong force gets smaller 
and smaller. Soon after this discovery [17>,18] . it was realized that this implies that at very 
high temperatures and/or very high densities the interactions between quarks grow smaller 
and smaller. These interactions could become small enough so that the quarks no longer 
bind together to form nucleons but instead roam freely through the system [|T5l . 
These free quarks have free color charges, just as the particles in a electromagnetic plasma 
have free electric charges. Hence, just as the charged particles in the traditional plasma 
screen the Coulomb interaction, free color charges screen the strong interaction (T3\. This 
screening washes out the long length scale interactions and the quark confinement is re- 
versed, what means that the quarks behave like free particles. At normal length and energy 
scales, free quarks are never observed and hence this prediction of free quarks is essentially 
a prediction of an entirely new form of matter being known as the Quark Gluon Plasma 

(QGP) m. 

Understanding the Quark Gluon Plasma could lead not only to new understandings of 
QCD, but also to new cognitions in cosmology and astrophysics. The conditions right after 
the big bang would have been ideal for the formation of a QGP, and there are also sugges- 
tions that color-charged matter might be important in neutron stars. 
In the recent years, a successful effort to solve the QCD equations numerically on a (space- 
time) lattice has brought deeper insight into the subject of phase transition(s) from hadronic 
to quark gluon matter ||T9l l20l |2T]| . So far it is not yet clear whether the transition is a true 
singular behavior of thermodynamic variables or just a rapid crossover. 
The only possibility of modern particle physics to produce a Quark Gluon Plasma is the 
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Figure 1.1: Result of a x^"fit to different particle ratios measured at ^/s=l30 GeV 



collision of heavy ions at ultra-relativistic energies. 

Heavy ion collisions 

The analysis of hadron yields measured in central heavy ion collisions from AGS up to 
RHIC energies has shown ^8311 that hadron multiplicities can be described very well with 
a hadro-chemical equilibrium approach governed by the chemical freeze-out temperature 
T, the fireball volume V and the baryo-chemical potential fiB- The RHIC result at -^5=130 
GeV is shown in Figure [TTTI for a certain set of particle ratios. 

The temporal evolution of a (central) nucleus-nucleus colhsion at ultra-relativistic energies 
is understood to proceed through the following stages [|23l : 

1. Liberation of quarks and gluons due to the high energy deposited in the overlap 
region of the two nuclei. 

2. Equilibration of quarks and gluons. 

3. Crossing of the phase boundary and hadronization. 

4. freeze out. 



at RHIC jni. 
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Reaching an equilibrium state is interesting, as it opens the door to a thermal description. A 
natural question though is how this equilibrium is achieved. Considerations about collision 
rates and timescale of the hadronic fireball expansion [[24| imply that at SPS, RHIC, and 
LHC the equilibrium cannot be established in the hadronic medium and that it is the phase 
transition which drives the particles densities and ensures chemical equilibrium [94J. 
Applying a statistical model which assumes an equilibrium and testing experimental data 
against model predictions is one way of testing reality against a thermally and chemically 
equilibrated fireball at the point of hadro-chemical freeze-out. The equilibrium behavior 
of thermodynamical observables can be evaluated by averaging over statistical ensembles. 
The equilibrium distribution is thus obtained by averaging over the complete accessible 
phase space. Furthermore, the ensemble which corresponds to thermodynamic equilibrium 
is the one for which the phase space density is uniform over the accessible phase space. In 
this sense filling the accessible phase space uniformly is both a necessary and a sufficient 
condition for equilibrium [231. 

e+e" collisions 

The discussion above showed that there is a good understanding of the system behavior 
in heavy ion collisions. An interesting question arising here is whether this statistical be- 
havior is a unique feature of high energy nucleus-nucleus collisions or whether it is also 
applicable in elementary collisions like, e.g., e+e^. Previous publications indicated that in- 
deed hadron production in e+e" collisions at 14-43 GeV [|6l|60l and 91 GeV [|6l|9T| can be 
described well within a Thermal Model provided that local quantum number conservation 
is properly implemented. The main results of these investigations are that the temperature 
values deduced are almost constant near T=160 MeV and that the volume increases with 
energy, while strangeness is under-saturated. These results were taken, together with the 
results for nucleus-nucleus collisions where a similar temperature is reached at high ener- 
gies, as evidence for the interpretation that the thermodynamical state is not reached via a 
dynamical equilibration among constituents. It is rather a generic fingerprint of hadroniza- 
tion [ | 24. ,92 J or a feature of the excited QCD vacuum L93i . 

In this thesis the latest multiplicity measurements of e+e~ collisions at ^/s=lO GeV, 29-35 
GeV, 91 GeV and 130-200 GeV summarized and published by the Particle Data Group 
(PDG) (5] are used. Since the aim is a precision calculation for a small system, a fully 
canonical form of the Thermal Model ll23l[T0l is applied, conserving baryon number, elec- 
tric charge, strangeness, charmness, and bottomness. In order to reach a precision compa- 
rable to that of the data from the LEP collider (a few percent) within the model, quantum 
statistics is included in the formulation. 

In chapter[2]the basics of statistical mechanics and thermodynamics are briefly introduced 
with focus on the equations which are relevant for this study. The following chapter [3] 
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1 Introduction 



discusses the understanding of an e^e~ collision, with a short overview over the exist- 
ing fragmentation models. The last part of this chapter introduces the Thermal Model, 
in the grand-canonical as well as in the canonical ensemble relevant for this work. The 
corresponding canonical partition function is derived in chapter |4l In contrary to earlier 
studies [161 13, the calculation of the partition function is carried out in quantum statistics 
instead of Boltzmann statistics. In chapter[5]the C-i-i- program, developed within this thesis 
to simulate the particle production in a collision, is explained and the main steps and ap- 
proximations are investigated. Section [6] demonstrates the influence of the canonical factor 
and shows the necessity for the canonical approach used in this work. In chapter |7] the 
fit procedure and the treatment of strong decays is explained. Furthermore the incidence 
of weak decays within a collision evolution is discussed and a significant influence to the 
particle yields is found. Finally all results are presented and discussed in chapter [H The 
appendix contains a fully quantum statistical derivation of the canonical partition function 
(for bosons and fermions) and a discussion of more technical aspects of the program de- 
velopment. 

The work presented in this thesis has been published in |[Tl|2l[3]0. 



' The results of this work and the resuhs pubhshed in the paper are not always identical. In such a case 
the values in the paper are the more recent ones. 



Units and conventions 

The common unit of length in hadron physics is 

1 fm = 1 Fermi = lO^^^m 

and for the time scale at relativistic velocities i> c 

Ifm 1 



lO^^^s 



c 2.99792458 
Energies are treated in units of 

1 MeV = 10^ • 1.602 • IQ-^^AsV = 1.602 • 10"^^! 

In all thermodynamic equations the Boltzmann constant is set to one 

kB = 1 

Therefore the temperature gets the dimension of energy 

1 MeV 1.60218 • lO-^^Nm ^ 

1 MeV = — = = 1.16044 • 10^°K 

ke 1.38066 • lO-^^Nm/K 

Also the speed of light is set to one 

c= 1 

and the Planck constant becomes 

h = 1.05457 • 10-^"Nms = 197.327 MeV fm 

1.60218 • 10-i3Nm W-^^m/c 



Chapter 2 

Statistical mechanics and 
thermodynamics 

In this chapter the basic properties of an ideal quantum gas are discussed. The equations 
presented in section [2H are used later to derive the partition function of a hadron gas. In 
section basic thermodynamical relations are introduced. 

2.1 Quantum statistics 

Particles of the same species are indistinguishable, for that reason it is not possible to 
specify the microstate of each particle. Instead a microstate of an ideal gas is specified by 
the occupation numbers Uk, the number of particles in each single particle energy Cfc. If the 
value of the occupation number is known, the total energy of a microstate \s) can be 
written as 

Es = ^nkek (2.1) 

k 

The set {n^} completely specifies one microstate of the system. The canonical partition 
function Z'^' of an ideal gas can be expressed in terms of the occupation numbers as 

s s 

= J2Y[e-'^'"''' (2.2) 

s k 

where (3 = 1/T0, V is the volume of the system and the occupation numbers satisfy the 
condition 

N = J2^k (2.3) 

k 

'in the rest frame of the system and fc^ = 1 



6 



2. 1 Quantum statistics 



7 



with the total particle number A^. 

According to the spin-statistics theorem all particles are classified into two groups. Par- 
ticles with zero or integer spin such as deuterons (^H^) or pions are bosons and follow 
Bose-Einstein statistics. Particles with half-integer spin such as electrons, protons, and 
neutrons are fermions and obey Fermi-Dirac statistics. The difference between fermions 
and bosons results in the allowed values for the occupation number n^. Fermions are re- 
stricted by 

rifc = or nfc = 1 (2.4) 

which is the Pauli exclusion principle. For noninteracting particles two identical fermions 
cannot be in the same single particle state. In contrast, the occupation numbers for 
identical bosons can take zero or any positive integer value 

nfc = 0,1,2,... (2.5) 

For fermions, the canonical partition function for the /c*^ energy level is 

1 

Z^{V, T)=J2 6"^"'-"'= = 1 + 6"^^'= (2.6) 

whereas for bosons 

(V, T) = J2 6-^"'=^'= = 1 + e-^'" + e-"^^'^ + ■ ■ ■ = ^ (g"^^^)"* (2.7) 

nfe=0 nfe=0 

This geometrical series is convergent for e~'^'^'= < 1, thus 

z-iy^T) = ^^^ (2.8) 

To get the average particle number at energy level k, one has to differentiate eq. (|2.6I) and 
(123 

These functions are called the Fermi-Dirac (-1-1) and Bose-Einstein (-1) distribution func- 
tions. 

For a system with more then one subsystems k, the partition function of the whole system 
is the product of all subsystem partition functions. This means that for fermions it can be 
written 

2f{V,T) = IIZ^{V,T) = + e-^^'=) (2.10) 



and for bosons 



Z^{V,T) = l[Z^{V,T) = 1113^ (2.11) 
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2 Statistical mechanics and thermodynamics 



In a similar way to the definition of the canonical partition function for the canonical en- 
semble, one can define a grand-canonical partition function for a grand-canonical ensem- 
ble. The only difference is a fugacity term 

z = e^^ (2.12) 

taking into account that the particle number in a grand-canonical ensemble is not fixed. 
The partition function of the grand-canonical ensemble is then given by a weighted sum of 
canonical partition functions with different number of particles N 

oo oo 
N=0 N=0 s 

For a system of bosons or fermions, it is often mathematically easier to treat the number of 
particles of the system as an intrinsic property of each state, |s) (specified by the set {uk}). 
Therefore the partition function can be written as 

Z'^^iz, V,T) = Y^ z^'^e-"''^ = n ^"'6-'^"'=^'= (2.14) 

s s k 

where A^^ is not fixed anymore. 

Now the discussion of the canonical ensemble (above) can be applied to derive the grand- 
canonical distribution functions 

Gc^ _ I gln-Zf^ _ 1 _ 1 
^'^'^ '~ 13 dpi ~ z-^e^'^ ± 1 ~ e/^(^*-'^) ± 1 ^ 

with — for bosons and -|- for fermions. 

2.1.1 Boltzmann statistics: The classical limit of Quantum statistics 

For fermions and bosons, the quantum statistical canonical partition functions can be writ- 
ten as 

lnZ^{V,T) = lnfnr3^)=-EMl-e-^-) (2.16a) 

\ k / k 

\nZ^{V,T) = In ( JJ(l + e-'^^'=) J = Zln(l-fe-^^'=) (2.16b) 

\ k / k 

where the results of the last section are used. The classical Boltzmann limit arises from the 
first term of the Taylor expansion of the logarithms 

ln{l±x)^^ + 0{x^) {x < 1) (2. 17) 



2.2 Statistical and thermodynamical relations 
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i.e., when it is possible to consider the exponential term as small compared to unity (e"'^"*^ <^ 
1). In this case the classical Maxwell-Boltzmann limit is reproduced for fermions and 
bosons: 

\nZ^i{V,T) = Y,e-^'' (2-18) 

k 

Furthermore it is also possible to apply the Boltzmann approximation directly to the distri- 
bution functions, introduced in (12.91) 

M = (2.19) 

where the plus sign refers to fermions, and the minus sign to bosons. The classical Boltz- 
mann approximation comes by neglecting the term ±1 in the denominator. Therefore the 
Boltzmann condition is again 

e^' > 1 (2.20) 

This is given for sufficiently low temperature or sufficiently high energy (mass) where for 
both, fermions and bosons, one obtains the same value 

M = e-^' (2.21) 

It follows again that the classical limit of the quantum distribution functions, for Fermi- 
Dirac as well as for Bose-Einstein, reduces to the Maxwell-Boltzmann distribution. 



2.2 Statistical and thermodynamical relations 

If the partition function of a system is known, it is possible to calculate all macroscopic 
properties of the system. For instance the average energy E is 

1 dZ d\nZ 

E = = (2.22) 

The free energy is defined as 

F = -T\nZ (2.23) 

and links statistical mechanics to thermodynamics. A pure thermodynamical definition of 
Fis 

F = E-TS (2.24) 
With the fundamental thermodynamic relation^ dE = TdS — pdV, one finds 

dF = -SdT - pdV (2.25) 



^The particle number conservation is not treated here 
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2 Statistical mechanics and thermodynamics 



Thus the entropy is 



and the pressure of the system is 



dF{V, T) 



P = 



dT 



dF{V, T) 



V 



dV 



(2.26) 



(2.27) 



In conclusion, the partition function Z is the key to all thermodynamical quantities. 



Chapter 3 

Theoretical description of e^e~ 
collisions 

The fragmentation process in an e+e^ collision is still an open issue, because a pure analyt- 
ical QCD approach is not able to access this area. This gave rise to many phenomenological 
models. After the discussion of the e~^e~ annihilation process some of these models are 
described. In this context the Thermal Model will be introduced. 

3.1 Characteristics of an e+e~ collision 

For a low center of mass (cm.) energy ^/s, the process e+e^ ff is dominated by a 
single virtual photon exchange. When the cm. energy of the resonance is reached 
(y^ ^ 91 GeV), exchange becomes the dominating process (see Figure [3121) . The most 
probable result of an e+e^ collision near the Z° resonance is multihadron production due 
to the large branching ratio of Z^ qq (^ 70%). Other channels, such as e+e~ — > e~^e~ 
(Bhabha scattering), e+e~ fJ'^l^^ and e+e~ r+r^, occur less frequently [l5l|. The 
further discussion will focus on the Z° decay to quarks, disregarding the leptonic decays, 
because this is the essential process needed in this work. 

The process e'^e" — > hadrons can be divided into four stages. In Figure [BTI these stages 
are shown: 

I. The formation of the Z^ in an e"^e^ collision, and its subsequent decay into qq are 
determined by the electroweak part of the Standard Model. 

II. As the quark and antiquark separate, gluons may be emitted, which may radiate fur- 
ther gluons or generate qq pairs. This process can be calculated in QCD perturbative 
theory, as long as the strong coupling, as, remains significantly smaller than one. 

III. The particle shower increases, with decreasing momentum exchange Q and increas- 
ing as- At some point as will reach the non-perturbative region {as ~ 1). Finally 
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3 Theoretical description of e+e collisions 
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Figure 3.1: The process e+e —> hadrons [1401 . 



the quarks and gluons will form colorless hadrons, but the transition from quarks 
to hadrons is not describable in perturbative theory, and one has to resort to various 
more or less ad hoc models to describe the process, usually called hadronization or 
fragmentation. 

IV. In order to compare theory with data one also has to count for decays of unstable 
hadrons in photons, leptons and stable hadrons. 



The electroweak stage 

In the first stage, the e+e^ pair annihilates into the resonance according to electroweak 
theory. This phenomenon may be accompanied by the emission of photons (initial-state 
radiation) prior to the annihilation. This electroweak correction reduces the cm. energy 
of the e"^e collision and therefore the total effective mass of the hadronic final state. 
Following its creation, the vector boson decays into a quark-antiquark pair (^ 70%). 
All these phenomena are described in the framework of the very successful electroweak 
model. 

The partial width for Z° — * qq is given by fi40il 

r,, = [(1 + 2r^^) {{gl.f + {9\f) - 6r^^(,^)^] (3.1) 



3.1 Characteristics of an e+e collision 
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Figure 3.2: Measurements of the different e+e^ cross sections at LEP. The mea- 
surements at lower energies are also shown for CESR, DORIS, PEP, 
PETRA and TRISTAN [193. 



with 

(5, = v^T^V; r^=^ (3.2) 

mz 

and with the vector and axial- vector currents 

gv = h - 2Q sin^ 9w (3.3) 
gA = h (3.4) 

For up-type quarks with Is=l/2 and Q=2/3 one gets 

9^ = ^-^sm''ew = 0.18 (3.5) 

9l = l (3.6) 

This gives r„=280 MeV. 

For down-type quarks with /3=-l/2 and Q=-l/3 one gets 

gt^ = ~ + ^sm^ew = -OM (3.7) 

9a = -I (3.8) 



14 
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Figure 3.3: Three jet event at the OPAL detector (LEP). 



This yields 1^=370 MeV. 

The total hadronic width from the three down-type quarks {d, s, b) and the two up-type 
quarks (u, c) is therefore rhadron=l-67 The relative hadronic branching ratios for the 

up-type quarks and down-type quarks are 



what will be an important point in the discussion of heavy particle production in e+e 
collisions, later in this work. 

The perturbative QCD stage 

In the second stage, the initial qq may radiate gluons according to the theory of Quantum 
Chromodynamics. The gluons may radiate other gluons or qq pairs, giving rise to a cascade 
process. This stage is responsible for the formation of hadronic jets as seen in Figure 13.31 
for a 3 jet event as it appears in the OPAL detector. The probability ratio for observing 
events with a certain number of jets can be written in terms of increasing powers of the 
strong coupling constant as [|33ll 



^uu /^hadmn ~ 17% 
r^^/Fhadron ~ 22% 



(3.9) 



2 jets : 3 jets : 4 jets = : 0{al) : 0{al) 



(3.10) 



'The top quark is not produced at -\/s=91 GeV. 



3.1 Characteristics of an e+e collision 
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Figure 3.4: The leading order Feynman diagrams of an e+e ^ Z ^ qq process. 



This means that at ^/s=9\ GeV the majority of the e"'"e~ events have a two-jet structure^. 
Due to the small value of a^, this stage can be described by perturbative QCD. Further- 
more the prediction of a cross-section a (either the total hadronic cross-section, or some 
differential cross-section) can be written in the general form 

a = aQ + Aas + Bal + Cal^ (3.11) 

To compute the leading coefficient, A, it is necessary to consider the Feynman diagrams 
shown in Figure [X4l The first two diagrams represent the radiation of a gluon from either 
quark or antiquark, leading to three-parton final states, and clearly yield cross-sections 
proportional to as- Although the square of the matrix element for the third loop diagram 
is O(a^), it leads to the same qq final state as the leading order diagram with no gluon 
emission, and the interference term between the two diagrams is 0{as) \56h 

The non perturbative QCD stage 

In the third stage, the partons fragment into colorless hadrons. Usually, this stage is called 
hadronization. This process cannot be described as a power expansion in the strong cou- 
pling constant, since > 1 (soft processes). 

Since the description of the final-state particles cannot be accessed analytically, a second 
approach is to use phenomenological models. Three of them are described later in this 
chapter. 

Particle decays 

The last stage in Fig. 13.11 represents the decay of unstable hadrons into experimentally 
observable particles (mostly pions). Many of these are well understood from low energy 
experiments, though in the case of heavy-flavored hadrons there are uncertainties since not 
all decay modes have been measured. 



^as at Vi=91 GeV is 0.12 ± 0.0031 ES). 
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3.2 Fragmentation models 

Quarks and gluons cannot be observed as free particles but are confined in color neutral 
hadrons. The process of hadron formation is not well understood, because it falls into the 
domain of non-perturbative QCD and it has not been possible to derive it from first princi- 
ples. Therefore, several phenomenological models aimed at describing quantitatively the 
fragmentation process have been developed. 

In quantitative terms, one has the picture of the color field lines between a quark and an- 
tiquark. As the quark and antiquark move apart their kinetic energy is transformed into 
increasing field energy, which can generate qq pairs, and turn into hadrons finally. There 
are basically three models: independent fragmentation, string fragmentation, and cluster 
fragmentation. The basic ideas of these models are summarized in the following subsec- 
tions. 



3.2.1 Independent fragmentation 

In the framework of independent fragmentation it is assumed that each parton hadronizes 
on its own in the following iterative scheme. A quark q fragments into a hadron qq' and 
a quark q'. The q' then proceeds to fragment in the same way, and the process is repeated 
until insufficient energy remains. The sharing of energy-momentum is governed by a frag- 
mentation function. Flavor and transverse momentum are conserved in each break-up. 
There are a number of difficulties with this scheme; for example, the result of independent 
fragmentation normally depends on the coordinate frame and is thus not Lorentz invariant. 
A fully Lorentz-invariant scheme has been proposed in [[42l| but it is hard to implement. 
Furthermore it does not conserve energy, momentum and flavor exactly (they are conserved 
locally but not global), so that these have to be patched up at the end. 



3.2.2 String fragmentation 

The quark and antiquark from the move apart in opposite directions. The color field 
between them forms a narrow flux tube called string. In the Lund model [|43l the string 
is idealized to an one-dimensional massless relativistic string with energy density k ~ 1 
GeV/fm. The potential rises linearly with the distance between the charges and at some 
point it becomes possible to form a new quark-antiquark pair from the field energy. 
In the classical picture, a qq pair with a mass or transverse momentum relative to the string 
is produced a certain distance apart, with the field energy between the quarks transforming 
into the transverse momentum. In quantum mechanics, they are produced at the same point 
(to conserve flavor locally) and then tunnel out. The tunneling probability for a common 
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transverse mass m]_ = + p^, of the qq pair is 11741 




(3.12) 



The tunneling also ensures that heavy quark and strange quark production is suppressed 
within the string. Heavy quarks in e'^e" annihilation thus either come from perturbative 
gluon splitting or are the primary quarks from the decay. The ratios of the different 
quark flavors given by the tunneling probability are u : d : s : c ^ 1 : 1 : 0.3 : 10~^^. 
Mesons are produced when a quark and an antiquark from adjacent string breaks are com- 
bined. The baryon production is usually modeled by the diquark [,45,1 or the popcorn [44| 
mechanism (see the references for details). 

The different break-ups of the string are assumed to be independent of each other. It is 
therefore possible to start at one end of the string, make a break-up and be left with a me- 
son and a shorter string. The procedure can be iterated until the available energy is used up 
(the termination requires some extra treatment that will not be discussed here). 
The fractions z of the available energy, taken by the hadrons, are distributed according to a 
probability distribution function f{z). The constraint that the result should be independent 
of the choice from which end to start leads to the Lund symmetric fragmentation function 

f{z) oc z-\l - zyexp{-bml/z) (3.13) 

Other functions are availablel^. The parameters a and b are chosen in order to fit the data. 
In the case of gluon bremsstrahlung, the situation is more complicated. 
Overall, the string model presents a more consistent and covariant picture than the indepen- 
dent fragmentation model and is generally in better agreement with experimental results. 
The Lund string model is implemented in the PYTHIA ll53l event generator (merged with 
the former JETSET B61I generator). The most recent version is PYTHIA 8.1 15411 . One 
drawback of PYTHIA is the large number of variable parameters, but it is considered to be 
the most faithful simulation currently available^. 



3.2.3 Cluster fragmentation 

The model of cluster fragmentation (4\] takes a very different approach compared to the 
string model. It is based on the observation that, at the end of the perturbative evolution of 
a parton shower, a given color and the corresponding anticolor generally occur on partons 
which lie quite close together in phase space. One can replace each gluon by a qq pair, and 
can then identify the color singlet combinations amongst all the quark and antiquark. Each 

^Such as the Peterson fragmentation function which is designed for the fragmentation of heavy 
quarks [i79l . 

'' ,http : / / home . thep . lu . se/~torb jorn/Pythia . htmT] 
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of these is considered to form a cluster, with flavor quantum numbers determined by the 
qq. The mass of a cluster is typically a few GeV. Each cluster is then allowed to decay into 
two mesons or a baryon-antibaryon pair, with the appropriated net flavor. 
This model typically has fewer arbitrary tunable parameters than the string model, but it 
gives quite a good account of the data. The most popular event generator which uses the 
cluster fragmentation approach is the HERWIG [55.1 event generatoi|f|. 

3.3 The Thermal Model of ultrarelativistic particle colli- 
sions 

The Thermal (Statistical) Model [l46l is a model of hadronization, aiming to reproduce the 
quantitative features of the hadronization process. This model has been recognized as a 
powerful approach to describe particle production yields in heavy-ion collisions [|23ir7TI . 
To be able to apply the framework of thermodynamics to a system of two colliding objects, 
the creation of a fireball in complete thermodynamical equilibrium at the collision point is 
assumed. The basic idea of the model is very simple and lies in two assumptions 1.46 J : 

1 . In the late stage of a high energy collision, the cluster (fireball) decays into hadrons 
in a purely statistical fashion. The decay happens at a critical temperature Tc of the 
fireball. 

2. All multihadronic states within the cluster compatible with its quantum numbers are 
equally likely. 

The analysis of particle yields measured in central heavy ion collisions from AGS up to 
RHIC energies has shown |[Ml[83[83[M[83[M[871IM|83l90l|7ll that hadron multi- 
plicities can be described very well with a hadro-chemical equilibrium approach which is 
governed by the chemical freeze-out temperature T, the baryo-chemical potential fiB, and 
the fireball volume V. 

In the limit of high temperature and/or large system size, the grand-canonical treatment is 
adequate. However, if the number of particles in a collision fireball is small, either due to 
low temperature (e.g. in nuclei-nuclei collisions at small energies of a few AGeV) or due 
to small system size (elementary particle collisions), then the canonical ensemble has to be 

used (a Ha. 

3.3.1 Grand-canonical ensemble 

In the grand-canonical ensemble the fireball is characterized by the temperature T (/3=1/T), 
its volume V and the chemical potentials /I which are assumed to be uniform over the whole 



• ^http : / /projects . hep forge . org/herwig/| 
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volume. These parameters determine the partition function 

Z^^{p,V,fI) = ^e-^^^^-'^^) (3.14) 

states 

In the "normal" grand-canonical approach the particle number N is not fixed, but conserved 
in the mean, by the chemical potential fi=iJ,^N. At relativistic energies the particle number 
is not fixed at all (because particles can be produced from the available energy). Instead 
of conserving the particle number, the chemical potential is used to conserve the quantum 
numbers of the system. Therefore /T is a combination of the different quantum number 
chemical potentials, e.g. /2=(yU7v, l^h^l^c^ f^s) for the baryon number N, strangeness 
S, isospin L^, charmness C and bottomness B. 

Note: It is crucial to understand, that the words grand-canonical and canonical in 
our context refer exclusively to the quantum number conservation laws, because the 
particle number is not conserved, neither in heavy ion, nor in e+e collisions. 

In the hadronic fireball of non-interacting particles and resonances, the partition function 
Z*^*^ is the product of the particle partition functions of all particles j (see chapter[2l) 

Z''^{(3,V,fl) = l[z]{(3,V,fi) (3.15) 
j 

with, 

zjiP,V,fI) = ^ r p^dp [l±exp(-/5(e,-x,-/i))]^^ (3.16) 
'^'^ Jo 

where ej is the relativistic energy of particle species j and gj={2Jj + 1) is the spin de- 
generacy factor. The vector Xj={Nj, Sj, I^j, Cj, Bj) contains the quantum numbers of the 
particle and jI={nN, /^5, ^ii^-iHc, (J-b) is the vector of chemical potentials related to the dif- 
ferent quantum numbers. The upper (lower) signs refer to fermions (bosons). 
The partition function (13.151) contains all information to obtain the particle multiplicity 
{nf^) of particle species j (see section [ITI) 

^ ' ^ 27T^J, exp[/?(e,-f, -/i)]±l ^ • ^ 

For hadron j with baryon number Nj, isospin I^j, strangeness Sj, charmness Cj and bot- 
tomness Bj, the particle specific chemical potential fij is 

/ij = Xj ■ fl = fiNNj + fiijsj + l^sSj + jJcCj + jj^sBj (3.18) 

The chemical potentials related to baryon number (/xat), isospin (fii^^), strangeness (/i^), 
charmness (fic) and bottomness (/i^) ensure the conservation (on average) of the respective 
quantum number: 

baryon number: V ^ rijNj = Z + N, (3.19) 
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strangeness: V^njS'j = 0, (3.20) 

■j 

isospin: ^^^^13^=/*°* (3.21) 

i 

charmness: V ^njCj = Q (3.22) 

i 

and bottemness: V ^ n,jBj = (3.23) 

j 

Where the (net) baryon number Z + N (Z = number of protons, = number of neutrons) 
and the total isospin /g"* of the system are input values which need to be specified accord- 
ing to the colliding particles (initial system) studie(i^. 

The temperature T, the volume V and the baryo-chemical potential are the three inde- 
pendent parameters of the model, while the chemical potentials /is, /x/3, /xc, are fixed 
by the conservation laws above0. T, V and are obtained from fits to experimental data. 
When fitting ratios of hadron yields instead of pure yields, the parameter V will cancel out. 



3.3.2 Canonical ensemble 

The grand-canonical ensemble is the simplest realization of a statistical approach and is 
suited for systems with a large number of particles. However, for small systems and for 
low energies, a canonical treatment is mandatory. It leads to a phase space reduction for 
particle production so-called "canonical suppression". 

The exact conservation of quantum numbers replaces the chemical potential of the grand- 
canonical ansatz by delta functions 

Z^{V,(3) = Y,e~^''S^^^^ (3.24) 

states 

where X = {N, S,Q,C, is the initial quantum number vector including all quantum 
numbers which are conserved and Xs contains the quantum numbers of the state The 
delta function forces the exact conservation of quantum numbers. 

The particle multiplicity in the canonical approach is suppressed, due to the fact, that all 
quantum numbers have to be conserved exactly (will be discussed in detail in chapter [6l). 
This canonical suppression is the main difference between the canonical and the grand- 
canonical ensemble at relativistic energies. A detailed derivation of the canonical partition 

''The degree of stopping of the coHiding particles, which is energy dependent and cannot be precisely 
determined experimentally, brings some uncertainty in the choice of initial baryon number and /g"*. 

^These conservation laws apply strictly only for quantities which are evaluated over the complete phase 
space. 

^In the canonical case it is often easier to use the electric charge Q instead of the isospin quantum number. 
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function is given in chapter El 

The canonical ensemble will find the grand-canonical results as the asymptotic value for 
large volume and/or large temperature (see chapter[6l). 



3.3.3 Strangeness suppression mechanism 

In the application of the Thermal Model to particle production in heavy-ion and particu- 
larly in elementary particle collisions it was found that canonical suppression alone is not 
sufficient to quantify the observed strange particle yields. Consequently, additional sup- 
pression mechanisms were proposed to account for deviations from experimental data. 
The suppression of strangeness has been parameterized by a factor, 7^, that is introduced 
to suppress hadrons composed of strange and/or anti-strange quarks 117211731 . 7^ is a value 
between and 1, where 1 means that all strange particles fill there whole phase space. In 
this description the particle partition function (eq. (13.161) in the case of the grand-canonical 
ensemble) of a particle composed of Ug strange quarks/antiquarks is modified in the grand- 
canonical as well as in the canonical ensemble by, 

4 ^7^ 4 (3.25) 

where 7^ is an additional parameter of the model and Ug. is the absolute strangeness of 
particle j. 



3.3.4 Thermal Model and e+e collisions 

The reach of an equilibrium is the reason, why the Thermal Model can be used to describe 
the hadron production of a high energy heavy ion collision. The question is whether this is 
also applicable to an e+e" collision, since it is a quite small system compared to heavy ion 
collisions. 

The result of previous investigations ^ l60l |9T1 was that the temperature values deduced 
are almost constant near T=160 MeV and that the volume increases with energy, while 
strangeness is under-saturated. These results were taken, together with the results for 
nucleus-nucleus collisions where a similar temperature is reached at high energies, as ev- 
idence for the interpretation that the thermodynamical state is not reached by dynamical 
equilibration among constituents but rather is a generic fingerprint of hadronization [|24ll92l 
or a feature of the excited QCD vacuum [93J. 

Alternatively, it was argued in [94J that the quark-hadron phase transition drives the equi- 
libration dynamically for nucleus-nucleus collisions. An apparent equilibration in e+e^ 
collisions is not easy to explain in this latter approach. 

Another difficulty in context of e+e^ collisions is the jet structure. The fragmentation 
models of section |3^ treat the particle production in e~^e~ collisions not separately in two 
jets. In these models the jets are combined by a string or cluster system, means they are 
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not independent of each other. This is possible, because the particle production is treated 
on parton levej^. 

If the Thermal Model is applied to an e~^e~ collision the assumption of the creation of a 
fireball in thermal equilibrium is needed. If only one fireball is assumed, it is hard to argue 
how a two jet structure can be produced. Therefore, the model usually works with the 
creation of two fireballs. N jet events caused by gluon radiation are neglected, which is 
reasonable (see section ISTI) . 

A further assumption is, that a jet keeps the original quantum numbers throughout its evo- 
lution. At this point a distinction is drawn between the so called uncorrelated and the 
correlated jet scheme: 

• In the uncorrelated jet scheme, both jets have to conserve the initial quantum num- 
bers. In an e+e" collision all relevant initial quantum numbers are zero. 

• A more general approach is to allow the sharing of quark pairs, means only both 
jets together conserve the quantum numbers of the initial system. This is called the 
correlated jet scheme. 

In this work both schemes, the correlated and the uncorrelated jet scheme, are treated. In 
the case of the correlated jet scheme it would make sense to transfer the quantum numbers 
of one quark originated by the decay, to one fireball. But then the quantum numbers of 
the quarks have to be conserved by the different fireballs separately. This will fail, because 
of the fractional quantum numbers of quarks like e.g. baryon number 1/3. The Thermal 
Model does not know any quark structure^, it knows only particle properties like mass 
and quantum numbers of the different hadrons and uses the Boltzmann factor to calculate 
their production possibilities. A fractional baryon number cannot be conserved through 
the production of hadrons (same for the electric charge). In addition to that the color is not 
conserved. This means the jets must interact at some stage of the hadronization process. 
The problem of fractional quantum numbers will be discussed in chapter[8l 



^If there is no treatment of separate jets in string or cluster models, how do they explain the jet structure in 
the final state? The jet structure is a result of the first quark pair arising from the Z'^ decay. For the radiation 
of hard gluons (gluons with much energy) small angles are preferred. The propagator factor is 

1 1 

(3.26) 



p2+p2 2E,Eg{l-cose. 



This means there is a preferred energy direction (= jet), the direction of the original quarks (quarks from 
the decay). This leads to a jet structure. 

"'One exception is the parameter 7^ which is discussed above. 



Chapter 4 



The canonical partition function at 
relativistic energies 

The exact treatment of quantum numbers in statistical mechanics has been well established 
ll63l [64l [65l [TOl for some time now. It is generally obtained by projecting the partition 
function onto the desired values of the conserved charge by using group theoretical meth- 
ods [|66l[67il . In this section these methods are introduced and it is shown, how one gets the 
partition function that accounts for exact conservation of quantum numbers. 
Later the derivation of a numerically manageable quantum statistical partition function 
is shown in detail. The partition function includes the conservation of baryon number, 
strangeness, charmness and bottomness. In section 14.2.11 the equations are extended to 
conserve also the electric charge. 

The results of this chapter are the equations used in the C++ program, developed within this 
thesis, to simulate the Thermal Model. All the results presented in chapter [8] are obtained 
by using this program. 

4.1 Canonical description of an internal symmetry for Abelian 
charges 

The description of a system in a canonical ensemble is based on a fixed number of parti- 
cles. However, in particle collisions at relativistic energies the number of particles is no 
longer conserved, as particles can be produced from the system energy. So the concept of 
the canonical distribution should be modified. 

The canonical distribution used here is related to the given irreducible representation of the 
symmetry group. Consider, for example, a gas consisting of particles all of which trans- 
form under the same representation U (g)". The transformation properties of the gas will be 
given by some reducible unitary representation U (g) which is a multiple tensor product of 
the representations U{g)°'. One can obtain states transforming under different irreducible 
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unitary representations of the symmetry group by decomposing the representation U {g) 
into irreducible ones. 

The usual way of treating the problem of quantum number conservation in statistical 
physics is by introducing the grand-canonical partition function 

Z{fik,V,p) =Tr[e-^(-^-^^'^^^'=)] (4.1) 

with the Hamiltonian H of the system, the charge operator, f3 the inverse temperature 
and the chemical potential associated with the conserved charge k, e.^. baryon number 
or strangeness. By denoting the states under the trace as \s) such that H\s) = Es\s) and 
Xk\s) = Xks\s) one can rewrite eq. (|4.1I) to 

Z(/ifc,l^,/5) = ^e-''^^ 6"'^^^^'=^^= (4.2) 

states 

The chemical potentials fik play the role of Lagrange multipliers which are fixed by the 
condition that the average value of the corresponding quantum number of a thermodynam- 
ical system is conserved and has the required value 

dfik 

This method, as is shown in chapter[6l is only adequate if the number of particles carrying 
the quantum number is asymptotically large and their fluctuation can be neglected. 
For simplicity only one quantum number, the baryon number A^, will be treated in the fol- 
lowing equations, in order to explain the transition to exact quantum number conservation. 
The generalization to more quantum numbers is straightforward. 
Then equation (14.21) becomes 

oo 

Z{f,^, V,(3) = Y, e-'^^^ e-^'^-^^ = ^ Z^{V, f3)X^ (4.4) 

states N=—oo 

where the fugacity Xn = exp(/5/iiv) is introduced and 

27v(^,/3) = 5]e-^^^- (4.5) 

SAT 

where sat is restricted to those states that carry an exact value of the conserved charge. 
Zj\f{V, (3) is the canonical partition function with respect to baryon number conservation. 
The goal is to calculate Z^iV, [3). 

Equation (14.51) can be rewritten with a projection operator 

Zr^iV, f3) = J2 e-"^- = E ^''""'P^ (4.6) 

Sjif states 
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where = is the projection operator on the states with exact value A^. For an Abelian 
symmetry, P^ is the 5-function, P/v = Sn,Ns 

ZN{V,f3) = J2(^~^'''^N,N. (4.7) 

states 

The conservation of additive quantum numbers like baryon number, strangeness, electric 
charge, charmness or bottomness is related to the invariance of the Hamiltonian under the 
Abelian U{1) internal symmetry group. In many applications it is important to general- 
ize the projection method to symmetries that are related to a non-Abelian Lie group. An 
example is the special unitary group SU (N) that plays an essential role in the theory of 
strong interactions. A generalization of the projection method would require to specify the 
projection operator. 

However, as all quantum numbers which are important for this work, namely baryon num- 
ber, strangeness, electric charge, charmness and bottomness are Abelian, the generalization 
to non-Abelian symmetries will not be discussed further, but the reader is referred to the 
original papers of Redlich and Turko lHH [6ll [62l| or a nice summery, given in [|23]| . 

4.2 Derivation of the partition function: 

Conserving baryon number, strangeness, charmness 
and bottomness 

In the case of the conservation of the baryon number, strangeness, charmness and bottom- 
ness the symmetry group is U{1)^, where each f/(l) corresponds to one of the conserved 
quantum numbers. Generalizing equation (14.71) to four exactly conserved quantum num- 
bers yield^l 

Z^,s,cAX) = Y.^~^'''^x,x. (4.8) 

states 

with the four dimensional vector X = {N, S, C, B). ^ is the usual Kronecker tensor, 

which forces the sum to be performed only over the states \s) whose quantum numbers Xg 
are equal to the initial set X. 

The integral representation of the delta-function above is 

-| 1*211 p2n i*2n /'2n 

'^^■^^ikL L 1 1 

'From now on the dependence of Z on 1/^ and (3 is not indicated anymore; instead the dependence on the 
quantum number vector X is pointed out. 
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Inserting (1431 ) in equation (|4~8]) yields (61 

1 /'27r /'27r /'27r nliv 

2n,s,cAX) = j^J / / / dVe^^^X^"'""^"^^^' (4.10) 

I^TTJ Jo Jo Jo Jo ,t,te, 

with = (07V5 05, 0C! 0b)- Each integral over one of the 0's represents the conservation 
of one quantum number. 

A state is specified by a set of occupation numbers {fij^k} for each phase space cell k and 
for each particle species j. This means Xs = J2j k where Xj = {Nj, Sj, Cj, Bj) is 

the quantum number vector associated to the j*^ particle. After summing over these states 
by using the results of section [2n eq. (14.101) becomes 

(4.11) 

with = number of bosons and Np = number of fermions. It is possible to rewrite this 
equation to 

Nb 



+ ^^^^ + e-'^"'^-'^'^) 

j=i k 

To get the mean number (multiplicity) of any hadronic particle one has to assign a fictitious 
fugacity Xj which multiplies the Boltzmann factor exp(— /Je^j — ixjcf)). This fugacity is set 
to one after differentiation 

, , d\nZ 



(4.13) 

A,=l 



hence 

1 ^ r f 

(n,) = / d^<P e^'^'^exp <^ V V ln(l 



'^n,s,c,b{X) (27 

Np 



(4.14) 
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where the + is for fermions and — for bosons. 

In the rest frame of the system, the four-vector (3 reduces to 

/3=(1/T, 0,0,0) 



(4.15) 



with the temperature T of the system. Furthermore, the sum over phase space cells in eq. 
(|4.12l) can be turned into an integration by changing to the continuum limit, where each 
quantum state refers to a classical phase space volume of (27r^)^ 



E 



9j 



dV 



(4.16) 



with dV = d^q d^p and the spin degeneracy gj. In a homogeneous system the volume 
integral Jy d^q gives a factor V and (|4.16l) is then 



E 



9,V 



d^p 



(4.17) 



(27r^)3 

The energy of a relativistic particle under the condition of a non-interacting system is 

= ^Jpl + (4.18) 
By using (I4.15L (14.171) and (I4.18L eq. (14.121) becomes 



Z 



N,S,C\B 



exp 



(27rfi)= 



(i^pln(l ± e 



-iXi(f)\H 



(4.19) 



For exp(— yp^ + m'j/T — zxj0) ^ 1 the Boltzmann approximation discussed in section 
I2.1.1l can be usecl^ 



ln(l ± e T 
Inserting (|4.20l) into (14.191) yields 

^0 



——ixjip 



Zn,s,c,b 



(2vr) 



(4.20) 



(4.21) 



3 



^For kaons the Boltzmann approximation causes an error of about 1%, but for pions the error can rise over 
10%. Therefore all particles with x ~ (0, 0, 0, 0) are calculated in correct quantum statistics (see eq. (I4.23l l). 
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with the particle partition function 



9^V 



n 2 I ™ 2 



and 



Zq = exp 



V 



{27rhy 



(fp e 



d^p In 1 - 



(4.22) 



(4.23) 



is the exponential of the sum over all pion-like states j, whose relevant quantum numbers 
are zero (states with xj = 0). With (14.131) and (14.211) one gets 



(rij) = 

whereas for all particles with Xj = 

{rijixj = 0)) = 



l2N,S,C,BiX — Xj 



z 



N,S,C,B 



(4.24) 



' a p- 



e T 



Assuming that ^ . z] <^ 1 for charm and bottom hadroni^ exp{^ 
can be expanded in power of z]^^^^ up to first order [[61 



(4.25) 



z] e-^v.''} 

In Ih J 



Jc/b ]c/b 



c/b 



^-l-^Jc/b"! 



ic/b ) jc/b 

Inserting (|4.26l) in (|4.21l) (for charm and bottom hadrons) one obtains 



Zn,s,c,b{X) 



(27r 



1 



3c 



= (27r)4 

1 Zq 



(27r)4 



ib & 



1 .1 ^0 



^"(271)4 



3c n ^ 



(4.26) 



(4.27) 



^The particle partition function zj is 0(10 ^) for charm particles and 0(10 for bottom particles. 
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with 

fi^) = J2^j^~''''^ (4.28) 
j 

where the index j runs over all charged hadrons except heavy flavored ones. The index jc 
runs over all charm hadrons and jb runs over all bottom hadrons. In this equation there are 
27r -periodic functions (the exponential functions) within 27r-integrals. The contributions of 
the charm and bottom particles are zero, except if the argument of the exponential function 
is zero. Therefore the integrals over (pc and 0^ can be replaced by Kronecker deltas 



jc ^ 



(4.29) 



Now, X and Xj are two-dimensional vectors having as components the baryon number and 
the strangeness, whilst the charmness C and bottomness B appear only in the Kronecker 
deltas. 

The integral representation of the partition function above is not convenient for numerical 
analysis as the integrand is a strongly oscillating function. However, it is possible to rewrite 
the integrals into sums over Bessel-functions [i7l[8l|9l. 
The integrals which have to be solved, are 

:Cn,s = d'^e'^^hxp jj] 4e-^^^^*^| (4.30) 

and after rewriting 



1 z-^' 


r 

d(j)N e'^'^- 




d(f)s e'^^^ i 


exp[ZA( 


' Qi{4>N-(t>s) _ 


exp [Zh ( 




exp[Zn( 


'g«(</'iV-3</)s) 



exp[Zx(e^'^^ + e-''^^)] 

_|_ g-«(</>]v-</'s)^j 
_^ g-«(</']V-2(/)s)^j 
_|_ g-«(</'iv-3(/>s)^j 



(4.31) 
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Quantum numbers Notation 



A^=0 


S=0 


(Zo) 


N=0 


S=l 


Zk 


N=l 


s=o 


Zp 


N=l 


s=-\ 


Za 


N=\ 


S=-2 


Ze 


N=l 


S=-3 


Zn 



Table 4.1: The possible quantum number combinations for particles (A^ = baryon 
number, S = strangeness). 



where ZK,p,A,E,n are sums over single particle partition functions having the same quantum 
numbers (see Table |4TI) . Note that all pion-like particles are already separated in Zq (see 
eq. (14.231) ). For the last three terms in eq. (|4.31l) the following relation of the modified 
Bessel-function can be used iTTSl 



exp 



J2 ^"^™(^) 



m=— oo 



and (14.311) becomes: 

oo oo oo 

n3=— oo 



By using 



one gets 



ni=— oo 



— / acpN e 

/TT 



?12=— oo n3=— oo 



1 

2^ 



27r 



COS( 



ni=— 00 

27r _ 
2^ 



n2=— 00 713=— 00 

i{N+ni+n2+n-i)(l>N 



exp[2Zp cos( 



'AT 



>s e 



i(S-ni-2n2-3n3)</)s 



exp [2Zk cos 05] 



(4.32) 



(4.33) 



(4.34) 



(4.35) 



4.2 Derivation of the partition function 



31 



Applying the following integral representation of the Bessel-function of order h 



1 f^"" 

Ih{x) = — / (i0exp(x cos 0) exp(— i/i^) (4.36) 
2vr Jq 



yields 



ni=— oo n2=~oo 713=— oo 

+ni+2n2+3n3 



(4.37) 



Now, the integral can be solved by summing over Bessel-functions. The last step is to 
replace the integrals in (14.291) by (14.371) to obtain 



jc 

+ Y i^N-N,,,s-sJc,oSB,B,, (4.38) 



3b & 



+ YY ^jc^jM-^Jc-N,,,S-S,,-S,^5c,C,jB,B,^+B,J 
jc n & ^ 

In the next section the charge conservation is included. The way used to derive a com- 
putable equation for Z{X) is very similar, therefore only the differences are shown. 

4.2.1 Model extension: Including conservation of electric charge 

The general canonical partition function is (eq. (14.191 ) for five quantum numbers) 

Zn,s,q,cAX) = ^ J d'^e'"^^ 



where X is now a five dimensional vector, X = {N, S, Q, C, B). To apply the Bessel 
relations (see (14.361) and (14.321) ) to the equation above, again Boltzmann approximation 
has to be used. But because it is now impossible to exclude the tt^ from the integrals (see 
the definition of Zq in (14.231) ) the Boltzmann approximation will give rise to a significant 
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deviation at least for the pion yield. 

The Boltzmann approximation is just the first term of the series of ln(l ± x)^^ (see eq. 
(|4.20l) ). To calculate it accurately (in quantum statistics) one has to use the whole series of 
the logarithm0 

oo ^ 

ln(l - x)-^ = TT (4-40) 

rh 

k=l 

yields 



Z 



N,S,Q,C,B 



(27r)5 

( oo 

exp \ ^ z]e-^^^^ + 5Z ^If^'"^^"^ 



(4.41) 



h k=2 



with 



where all particles with x = are excluded within Zq, which is defined in (|4.23l) . but here 
just for 7r°-like particles. The index j runs over all charged particles and b runs only over 
the charged light boson^. 

Applying (14.131) to (14.411) gives for bosons, e.g. the n"^ 



, , d\nZ 



A„±=l 



z ± ■ 



Zn,s,q,c,b{X) Zn^s,q,c,b{X) 

fc=i 2^n,s,q,c,b{X) 
whereas for fermions and heavy bosons eq. (14.241) is still valid. 

After the approximation for charm and bottom particles (see eq. (14.261) ) the following 



"^Because of the series used (eq. ( 14.401 )) the quantum statistical calculation is only applied for particles 
which follow Bose-Einstein statistics. The extension of the equations to fermions is straightforward (see 
appendix IHll. However, the deviation for fermions caused by the Boltzmann approximation is negligible 
because of the high mass. 

''Charged means particles with x ^ 0, not necessarily connected to the electric charge Q. 
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integrals have to be solved 

'^N,S,Q '■ 



After rewriting 



1 



(27r)3 
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2tt 



oo 



b k=l 
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2lT 



2n 
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1 

2^ 
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i(4>N-'Ps+<f>Q) 



+ e 
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_|_ g-*(0B-2</>s)-)j 
i{<t>N-'i<l>S-4>Q) _|_ Q-i{<t>N~'i<f>S 
i{<j}N-S4>S-<l}Q) _|_ p-i(0JV-30s 



-<?^q))] 



exp[ZnT(,e 

oo 

exp ^Z^±(e^'=<^Q +e 



+ e" 



exp 



exp 



k=2 

oo 



/i=2 

oo 



i=2 



(4.44) 



(4.45) 



where the Zs are given in Table 14 .21 Using the Bessel relations (14.321) and (14.361) gives 

10 oo 

/„,1 (2Zp)/„2 (2Za¥ )/n,3 (22'a++ )/n4 (2.^_ft:± ) 



X, 



N,S,Q 



n E 

j=l Uj =— oo 
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Quantum numbers Notation 
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Qr=0 


(^o) 


AT n 

iV =u 


C 1 


A 

v=o 




A'=l 


O A 

D=0 


Q=0 




AT A 

iV =u 
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r\ 1 

(^=1 




AT 1 

A'=l 


6=0 




ry 

Zjp 


iV=l 


<s=o 


Q=-l 


ry 


AT" 1 

iV=l 


O A 


Q=2 


ry 

^A++ 


AT" A 
iV =U 


O 1 

o=l 


1 




N=l 


S=-l 


g=o 




N^l 


S^-l 


g=i 




N=l 


S=-l 


g=-i 




N=l 


S=-2 


g=o 




N=l 


S=-2 


g=-i 




N=l 




g=-i 





Table 4.2: The possible quantum number combinations {N = baryon number, 5" = 
strangeness, Q = electric charge). 



In5{'^Z\)In6{2Z-£+)In7{2ZY.-)Insi'2Zso)In9{2Zs^^)Inloi2ZQT] 



oo oo 



n E 

_k=2 nfc=— oo. 



oo oo 



oo oo 



n E 

J=2 n(=— oo. 



/„,(2Zj^±) 



n E 

Ji=2 nh=—oo_ 
I-N-nl-n2-n3-n5-n6~n7-n8-n9-nlo{'2Zn) 
I-S-n4+n5+n6+n7+2n8+2n9+3nl0-J2 hnu-T, '"i (^-^K") 

-^-Q-nl+n2-2n3-n4-n6+n7+n9+nl0-X; kuk-'E^'^l (^-^tt*) (4.46) 

and the final partition function for the conservation of five quantum numbers is 



2n,s,q,c,b{X) fti Zo\^N,s,Q ^c,oSb,o 

+ 5Z ^jc^^-^jc,S-Sj,,Q-QjJc,CjjB,Bj, 



3c 
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36 & 



+ X! ^jd^^-^jc-Nj,,S-Sj,-S^^,Q-Q^,-Qj^^Sc,CjjB,Bj^+Bj^ j 

jc ib & ^ 



(4.47) 
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4.3 Alternative way to solve the canonical partition func- 
tion 

There is another possibility to calculate the canonical partition function (14.411) numerically. 
This way of calculation was first published in lfT4ll . but only for Boltzmann statistics, and 
with the conservation of three quantum numbers. Here the quantum statistical case with 
the conservation of five quantum numbers will be discussed. 

The calculation of [14] is implemented in THERMUS [98l, a Thermal Model simulator 

which is used later for a comparison to the results of this work (see chapter[8l). 

The canonical partition function with the conservation of five quantum numbers is (eq. 

dSB) 




(4.48) 



k=2 



with 



k{2nhy 



= TT^s I d'P e--^' (4.49) 



Remember that j runs over all charged hadrons, whereas b runs over the light charged 
bosons (also in this section the full quantum statistical calculation is treated for bosons 
only). Zq contains the contributions of all particles with xj = (see eq. (|4.23l) ). 
First the approximation (14.261) is applied to get rid of the integrals over (pc and (pB- Next, 
the exponent of the exponential function can be separated into mesons and baryons 



00 00 



E^^-^+EE'^^""'*- E E4- 



h I. — 9 mesons ju — 1 



baryons 



(4.50) 



baryons 



00 



+ ^ ^ 2z^ cos{kSj(j)s + kQ. 



mesons fc=l 



where the sums over mesons and baryons extend only over the particles (i.e. not the an- 
tiparticles). Defining 



uj= zje^(^^'^^+^^'^«) = \uj\e''''^'^ (4.51) 

baryons 
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the two sums over baryons can be written as 2|ci;| cos(07v + argcj). Substituting these 
expressions into the three remaining integrals of (I4.48|) gives 

■1 nlTZ nliT pliv 

^N,s,Q = TTTTT / d(j)^ e^^<^- / d(l)s e'^'^^ I d<j)o e 



(2vrj ,0 

exp(2|c^| cos(0iv + argu;)) (4.52) 

{oo 
2z^^ cos{kSj(l)s + kQj 
mesons j k=l 



Changing variables to </)'^ = </)jv + arg u gives 

l^^j JQ Jo 



and so 



exp(2|^| cos(0'^)) (■4 53) 

{oo 
mesons j 



(27r 

f.27r 

27r 



JO 



1 Z"^'' / 

— y (i0'^exp(2|u;| cos0'^) exp(2A^0'^) (4.54) 



exp 

I. mesons j k=l 

Using the integral representation of the Bessel-function (14.361) 



yields 



Ih{x) = — / (i0exp(x cos 0) exp(— z/i0) 
27r 







exp <{ 5Z 5Z ^4 coi^{kSj<j)s + fcQi</'Q) [ /-iv(2|cj|) 

^ mesons j k=l 



(4.55) 
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Finally, using the fact that the partition function is real (in fact, by symmetry, the imaginary 
part integrates to zero) 



^N,s,Q = -nz-TK 1 d(j)s I d(j)Q cos{S(j)s + Q(t>Q - Naigu) 
(27r)^ Jq Jq 

c oo . (4-56) 

exp J 5^ 5Z cos(fcS'j</)5 + fcQj^g) > /-7v(2|cj|) 



, mesons j fc=l 



Now a double integration is left which can be performed numerically. The complete parti- 
tion function with the conservation of five quantum numbers is again (|4.47l) 



J6 & 



jc n ^ ^ 

(4.57) 

A comparison of the calculation in section 14.2.11 to the calculation shown in this section 
is given in appendix lA.ll Both methods agree very well. Nevertheless in this thesis the 
calculation explained in section 14.2. II is employed and all computations later on are done 
this way. 



4.4 The final partition function 

Up to now only the partition function and the multiplicities within one jet are calculated. 
As already mentioned an e+e~ collision, emerges into two jets. Therefore the partition 
function of the whole system is 



Z = Z{X)Z{-X) 



(4.58) 
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Deriving this function gives the multiplicities: 



d\nZ 



1 dZ 



A,=l 



ZdXj 

l_ f dZ{X) 
Z\^ 
1 



Z{-X) + Z{X) 



dZ{-X) 
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'i V z{x) ^ z{-x) ) 



Z{X-kxj) _|_ Z{-X-kxj) 



Z(X) 



Zi-X) 



(for fermions) 
(for bosons) 



where the product rule and dZ{X) /dXj = ZjZ{X — xj) was used. 

In the case of a neutral system, like e+e~, all initial quantum numbers are zero (X = 0) 

and therefore the equation above becomes 



(n,(X = 0)) 



2;z, 



1 Zj-Xj) 

i 2(0) 

oo _kZi~kXj) 



1 Z{0) 



(for fermions) 
(for bosons) 



If sharing of quark pairs between the two jets is allowed (correlated jet scheme), the vector 
X is not equal to the zero vector, and one has to use the general form of the equations. 
Furthermore in the correlated jet scheme the multiplicity of particle j is the sum over the 
production in the different jets 



{nf') = Y.R,{n,iX,)) 



(4.59) 



Rg is the fraction of the specific jet (^ 17% for u-type quarks and ^ 22% for d-type 
quarks, see equation (|3.9I) ). The vector of quantum numbers Xg for a u-jet is Xu = 
(1/3, 0, 2/3, 0, 0), whereas for a 6-jet it is X^ = (-1/3, 0, 1/3, 0, 1). 



Chapter 5 

The C++ program 



The equations discussed in the last chapter are implemented in a class based C++ program 
developed within this thesis in order to calculate the hadron production in particle colli- 
sions. A detailed explanation of the program code can be found in a Bachelor thesis 
developed parallel to this Diploma thesis. Here, just a short overview is given. 

5.1 Structure of the program code 

The base of the program is a table storing all relevant particle properties i.e., mass, degen- 
eracy, resonance width, quantum numbers etc. The table includes all available particles up 
to a mass of 3 GeV and all important decay channels for each of them using the most recent 
PDG publications [5j. Each particle is a member of the class SParticle. See for instance 
the A+: 

Listing 5.1: A"*" particle as a member of the class SParticle 



case 233: { strcpy (name ," Delta +"); g = 4; m=1232; 



width = 120; ispin=0.5; 


charge 


=1; threshold 


=1073; 


baryon = l; strange=0; 


charm = 


; bottom=0; 




strcpy (channel[l] ," 


Neutron 


"); pchan[l] 


=0.333; 


strcpy (channel [2] ," 


Proton " 


); pchan[2] 


=0.667; 


strcpy (channel [3] ," 


Pi +"); 


pchan [3 ] 


=0.333; 


strcpy (channel [4] ," 


Pi 0"); 


pchan [4] 


= 0.667; break;} 



The A+ particle has a degeneracy of 4, a mass of 1232 MeV, a resonance width of 120 
MeV etc. The four most important decay channels are the decays to neutrons, protons, 
7r+ and 7r°. The corresponding branching ratios are listed in the array pchaii[]. The mass 
spectra of all particles collected in the particle table is shown in Figure ITTTl 
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Initial 
parameters: 
T, V and 7, 



Particle 
initialization 



Calculation of 



and lyiZ^j 



Calculation 

of Z{X) 



Build LUTl: 
Containing the 
Bessel-functions 




DecayMachine 



no 



Write to file 



Figure 5.1: Flow chart of the program used to calculate the particle multiplicities. 
For details see text. 



5.1 Structure of the program code 
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In the following the program behavior will be explained, going through the sequences 
which have to be passed in a normal run of the program as presented in Figure lSTTl The first 
step is the initialization of the particles using the particle table described above. Moreover 
the informations stored in the particle table can be used to calculate the arguments of 
the Bessel-functions (see Table |4T| and Table 1431) . Therefore the program considers all 
particles, collects their necessary properties and calculates the particle partition function of 
each of them via 



(5.1) 



After this procedure these functions are added up for particles with the same quantum 
number 

~ (5.2) 



Now the arguments of the Bessel-functions are known. During the calculation there are 
(9(10^) calls of the Bessel-function, but it is possible to reduce the number of calls to 
(9(10^) by using a lookup table (LUTl). The Bessel-functions are calculated by employing 
the routines published by "numerical recipes in C-i-i-"[.l 1,|. 

In the next step, the partition functions Z{X) and Z{X — kxj) are calculated in a loop 
over all particles by using equation (14.471) . The function is called ParticleCalc() and needs 
an array as argument, containing the quantum numbers X or X — kx. Now all necessary 
values are available to calculate the multiplicities 



.Z{X — kxj 



(5.3) 



where k is fixed to 1 for all fermions and heavy bosons (see chapter H]). 

To sum over Bessel-functions (from now on called Bessel blocks) in eq. (14.471) like for 

instance 



10 

n 

j=l rij 



oo 



Indi'^ZjC) In6{2Z-s+) Inli'^Zj^-) Insi'^Z^o) In9{2Zs^) Inloi^2Zfi^^] 



n E 

.k=2 nfe=— oo. 



n 



oo 



.h=2 nfi=—oo. 



oo oo 

n E 

.1=2 n;=— oo. 



(5.4) 



-^-Af-nl-n2-n3-n5-n6-n7-n8~n9-nlo(2.^n) 

-^-S'-n4+n5+n6+n7+2n8+2n9+3nl0-X; ^"h-E '"i (^^-f^'") 
-^-Q-nl+n2-2n3-n4-ji6+n7+n9+nl0-E knk~Y, (2.^7r± ) 

is very time consuming. During one run (9(1000) calls of this Bessel blocks are needed, but 
since the Bessel block depends on the quantum numbers only, it is the same for all particles 
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with identical quantum numbers. Therefore it is very helpful to build a second lookup table 
(LUT2) containing the results of the previous calculations of the Bessel blocks. This table 
is build during the runtime. If a Bessel block is needed, it should be checked whether the 
block is already calculated. The results are stored in a three-dimensional array with the 
quantum numbers as indices. Therefore each search in the LUT can be done with 0{1). 
This extension reduces the number of calls from (9(1000) to (9(10). 
After the calculation of all multiplicities, the decay of the heavy particles should be taken 
into account. The corresponding routine in the program is called Decay Machine(). It uses 
the branching ratios provided by the particle table and calculates the contribution of the 
heavier particles to the lighter ones. Only stable particles do not give decay contributions. 
The property "stable" is a definition depending on the measured particle yields. With the 
used measured data ^ all particles with a decay length cr > 10 cm are stable (see chapter 



In Listings |S!2I the first part of the output of the correlated calculation is shown: 

Listing 5.2: program output (only the first part) of the correlated jet scheme 
root(auser : / $ ./therm -cor 157 32 0.8 test.txt 

Calculating... T: 157 * V: 32 * gams: 0.8 ... 
*** Total nr. of particles: 474 *** 

We are in a correlated jet scheme 
weak decays : 1 

kpion : 5 
kkaon : 3 
kckaon : 3 



X[0]: 0.3333 

X[l]: 

X[2]: 0.6666 

X[3]: 

X[4]: 



m. 



Please wait and count the stars: 
round : 2 



Z(X) 



= 43.59470046778414 



5.2 Numerical and mathematical approximations 
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The routine needs some parameters, namely the scheme (-cor = correlated jet scheme), the 
model parameters T, V and 7^ and an output file for the results. The terminal output con- 
tains the number of particles (474) and the state of the weak decays (see section |7]). The 
three variables kpion, kkaon and kckaon give the indices k, h and / of the products in eq. 
(|5.4I) and will be explained in the next section. 

Although the program runs in the correlated jet scheme, the output in Listings 15.21 is re- 
stricted to the first jet, an w-jet, defined by the quantum numbers A^=l/3 and Q=2/3. The 
anti-M jet will be calculated in another loop with inverted quantum numbers. The result for 
the partition function Z{X) is also showri]. 

The numerical implementation of the equations was a big challenge, as the possibility of 
a fit procedure which is crucial for a data analysis depends strongly on the runtime of the 
program. With the implementation of the two lookup tables and a precise analysis of the 
necessary terms in the calculation process (see next section) it was possible to reduce the 
runtime to an acceptable value. In the correlated jet scheme the program takes 13 min- 
utes (on an Athlon64 2600MHz, LI 2x 128 KB, L2 2x 1024 KB, FSB 1066MHz) for one 
calculation. The other schemes are much faster. 



5.2 Numerical and mathematical approximations 



Equation (15.41) shows some problematic points for numerical implementations, e.g. terms 
like 



00 00 

n E 

,k=2 ?ife=— 00. 



00 00 

n E 

./i=2 nfi=—oo_ 



00 00 

n E 

1=2 n;=— 00. 



/„,(24.±) (5.5) 



This terms ensure the correct quantum statistical calculation for light bosons. But the math- 
ematical calculation needs upper limits for k, h and /. 

The contribution of the Bessel-function depends on the argument Z,j. and this argument 
decreases very fast with increasing k (h, I), because the index is included in the argument 
of the exponential function in (|5.1I) and is so exponentially suppressed. Therefore just a 
few terms labeled by the index k (h, I) are expected to be needed, but this depends on the 
model parameters. This is shown in Figure [5]2] for a certain parameter set (T=157 MeV, 
V=32 fm^ and 7^=0.8). It can be seen that for pions not more then five terms of the series 
of the logarithm (|4.40l) are needed to get the correct quantum statistical particle yield and 
even less terms for kaons. It is also visible that the deviation depends on the mass of the 
particles and decreases quite fast with increasing mass. For protons the deviation is much 
below 0.1%, what gives the reason to use Boltzmann approximation for fermions. 



' Zq is not included in the calculation, since it will cancel out later. 
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Index k (h, I) 

Figure 5.2: Deviation of the initial 7r+ and yields from the quantum statistical 
value with increasing index k (h, I). 
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Figure 5.3: Bessel-function for different particle species. 



5.3 Fitting with ROOT MINUIT 
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The next point deals with the infinite sums of eq. (15.41) like 



oo 



J2 ^n{2Z,) 



(5.6) 



n=— oo 



because also here a finite index n is needed. 

The biggest terni of the sum is the term with index 0. This is shown in Figure 15.31 for 
different Z^. With increasing or decreasing index n the value decreases quite fast. All 
Bessel-functions Ix{y), in the parameter range of this work, are not much larger than 1, 
which means that the smallest contribution in the product of (15.41) gives the upper limit to 
the dimension of the whole Bessel block. This can be used to set limits for the indices of 
the sums in (|5.4I) . 



To perform a fit procedure, the program is linked with ROOT MINUIT [HI. The MINUIT 
package acts on a multiparameter object function which will be minimized. This function is 
called the FCN function for which MINUIT defines the pure abstract base class FCNBase 
as interface. The FCN function must be implemented in a derived class from FCNBase. 
The value of the FCN function will in general depend on the parameters of the system. In 
our case these are the temperature T, the volume V and 7^. In the case of a fit procedure it 
is common to use a chi-square or maximum likelihood function to minimize, but in general 
this is a choice of the user. Here a chi-square function is employed: 



With this implementation, MINUIT is requested to minimize the FCN with respect to the 
parameters, that is, find those values of the coefficients which give the lowest value of 
[I103L This function will be discussed further in chapter|71 



5.3 Fitting with ROOT MINUIT 
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(5.7) 



Chapter 6 

About the necessity for the canonical 
approach in e+e~ collisions 

In chapter |4] the canonical description of the conservation laws in e"*"e~ collisions was 
used although in heavy ion collision a grand-canonical ansatz seems to be acceptable. The 
canonical description is valid in the whole parameter range whereas the grand-canonical 
formalism is the asymptotic realization of the exact canonical approach. However, the 
increase in complexity from the grand-canonical to the canonical ansatz is large (seen in 
chapter!!]) and rises the question whether it is necessary to use the canonical formalism, or 
whether it is allowed to perform the calculations within the much simpler grand-canonical 
framework. 



6.1 Grand-canonical conditions 

The particle number fluctuation in a grand-canonical ensemble is [fTOll 

with ^ oo (6.1) 

Hence for large the result of a grand-canonical ensemble is equivalent to a system with 
constant N, i.e. a canonical ensemble. Therefore canonical and grand-canonical ensemble 
give the same results at the limit of large A^. 

The system treated here has the temperature T and the volume V as free parameters. By 
relating T and V to the number of particles it can be investigated, at which limits for T 
and V the grand-canonical treatment is valid. 

The temperature T and the volume V are included in the particle partition function 



(N) Jn 



46 



6.2 The canonical factor 
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(particle phase space). Moreover Zj is proportional to the multiplicity 

~ ^1 ~ {rxp{-i/T} 

As one can see higher numbers in T or F increase the particle partition functions and 
therewith the multiplicity of the particles (A^). Therefore the following limits are equal: 

N ^ oo 

V ^oo (6.3) 

T ^ oo 

This is treated in detail in the next section. 



6.2 The canonical factor 

The grand-canonical approach gives the following equation for the multiplicity of particle 
species j fTTI 

oo 


with spin degeneracy gj, momentum p, total energy and chemical potential /ijQ. In the 
e+e~ system all chemical potentials are zero, and eq. (|6.4I) becomes 

oo oo 







where the Boltzmann approximation was used in the last step. The right hand side of 
this equation is identical to the canonical particle partition function of chapter HI equation 
(14:221) . 

The canonical approach gives for the multiplicities (see eq. (|4.24l) i^ 



oo 



(n^) ~ ^^^^ M I d^r) e-^'^ (6 6) 







'/ij is a combination of different chemical potentials, depending on the conservation laws treated. E.g. 
Mj — t^sBj — iisSj — f^hlj where the quantities Bj, Sj and /| are the baryon number, strangeness and 
three-component of the isospin quantum numbers. 

^For simplicity everything is treated in pure Boltzmann statistics. 



48 6 About the necessity for the canonical approach in e+e collisions 



o 

N n 



0.8 



0.6- 



0.4 



0.2 



.iiii.iiiiii. 



1 




10' 






10-= 






10' 






lO" 












10 = 






10-' 






10-= 






10-= 






lo-'- 






10-" 








1 1 1 1 





012340123401234 

baryon number strangeness electric charge 



1 
charmness 



1 
bottomness 



Figure 6.1: Development of the partition function with increasing baryon number, 
strangeness, electric charge, charmness and bottomness. 



Therefore the only difference between the grand-canonical and the canonical description is 
the canonical factor Z'"(X — Xj)/ Z'^(X) (in neutral systems). This "correction" factor 
depends on the thermal parameters of the system and the quantum numbers of the particles 
only (i.e. the correction for the A+ and the proton are the same). It is also obvious, that 
for neutral particles (xj = 0) the canonical factor equals unity and eq. (16.61) is exactly 
the grand-canonical result. The canonical factor is a typical feature of the canonical ap- 
proach due to the requirement of exact conservation of the initial set of quantum numbers. 
This factor suppresses or enhances the particle production according to the vicinity of their 
quantum numbers to the initial X vector. 

The behavior of Z'^{X) as a function of the baryon number, strangeness, electric charge, 
charmness and bottomness for suitable T, V and 7^ values is shown in Fig. 16.11 For in- 
stance, it is evident that the baryo-canonical factor Z^{N, 0, 0, 0, 0)/Z'"(0, 0, 0, 0, 0) con- 
nected with an initially neutral system plays a major role in determining the baryon multi- 
plicities. 

The ultimate physical reason of charged particle (xj 7^ 0) suppression with respect to neu- 
tral ones {Xj = 0) in a completely neutral system {X = 0), is that for every charged particle 
created, there has to be an anticharged particle, so that the conservation law is fulfilled. In 
a finite system this pair creation mechanism is the more unlikely the more massive is the 
lightest particle necessary to compensate the first particles quantum numbers. For instance. 
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neutral system 
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Figure 6.2: Volume and temperature dependence of the canonical factor for a 
neutron (antineutron) Z'^{1, 0, 0, 0, 0)/Z*^(0, 0, 0, 0, 0). It starts with 
T=150 MeV and goes - in 10 MeV steps - up to 200 MeV. The blue 
shaded area is the volume region investigated in this thesis. 



once a baryon is created, at least one anti-baryon must be generated, which is rather un- 
likely since its mass is much greater than the temperature and the total energy is finite. 
On the other hand, if a non-strange charged meson is generated, just a pion is needed to 
balance the total electric charge. Its creation is clearly less unlikely with respect to the 
creation of a baryon as the energy to be spent is lower. This argument illustrates why the 
dependence of Z'^(X) on the electric charge is much weaker than that on baryon number 
and strangeness (see Fig. 16.11) [80] . 

These canonical suppression effects do not occur in a grand-canonical framework. Equa- 
tion (16.51) shows that in a completely neutral system, all chemical potentials are zero and 
consequently charged particles do not undergo any suppression with respect to neutral ones. 
With increasing temperature the probability of the production of additional particles to bal- 
ance the net quantum number is increasing, too. Thus the canonical factor should converge 
to unity. This is shown in Figure 16. 2[ One can also see the conversion against 1 with 
increasing volume. This is already mentioned in section 16. 1[ with a slightly different ar- 
gument, but both explanations have the same origin, the increasing number of particles. 
Therefore the general behavior of the canonical factor at the limits found in section l6n is 
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Figure 6.3: Canonical factor of p, p, 7r+ and tt" in a proton-proton system. 

= 1 (in e+e^coUisions) 

and the system transits from the canonical to the grand-canonical ensemble. 
Figure [6^ shows an example for a baryon without strangeness or electric charge (neutron). 
It can be seen that the canonical factor contributes a non-negligible effect to the particle 
yield at the volume and temperature region discussed in this thesis (blue shaded area, see 
chapter[8]). This justifies the canonical treatment in chapterlH 

A detailed look at special particles in a neutral system, shows that there is no difference 
between the suppression of a particle and the corresponding antiparticle. The production 
of a baryon needs the additional production of an antibaryon, and the production of an 
antibaryon needs a baryon, so one finds 

Z^{N, S, Q, C, B) = Z^{-N, -S, -Q, -C, -B) (6.8) 

Therefore we have the same yields for particles and antiparticles in a neutral system. 

6.2.1 Charged systems 

In a "charged" system like for example the collision of two protons with the initial quan- 
tum number vector X = (2, 0, 2, 0, 0), it is not possible to use the direct analogy between 
canonical and grand-canonical systems as done in the neutral case in eqs. (16.51) and (16.61) . 
because the chemical potentials are not equal to zero anymore. Nevertheless the suppres- 
sion effect of the canonical factor can be investigated. But this effect changed now dra- 
matically for some particles. Figure [63] shows an enhancement for protons, instead of a 
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220 



volume[fm ] 

Figure 6.4: Canonical factors for K° and K° in the proton-proton system. The K° is 
more suppressed than the K°, as it is possible to use baryons to balance 
the strangeness for the K° (see text). 



suppression. The explanation is the following: The canonical factors are still intended to 
conserve the initial quantum numbers. To conserve a baryon number of two, more baryons 
and less antibaryons are needed. The same effect can be seen at the pion yields. The tt"*" is 
enhanced, due to the positive electric charge in the initial system. But pions are produced 
easily, therefore the effect is much smaller than for protons. 

It is also interesting to look at the canonical factors of K° and K'^ in a proton-proton sys- 
tem. The initial system contains no strangeness, therefore one would expect for both kaon 
canonical factors the same behavior as in a neutral system. The canonical factors are plot- 
ted in Fig. 16.41 There is a suppression for both - kaon and antikaon - but the suppression 
for the antikaon is larger. This can be explained by looking at the details of quantum num- 
ber conservation. To balance a quantum number, it is convenient to look for the lightest 
particle which can do the job, because this particle has the biggest probability. But in a 
pp system it has to be taken into account, that baryons are enhanced and anti-baryons are 
suppressed. To balance the strangeness S = 1 arising by the production of a K° it is still 
the easiest way to produce a K°. But its also possible to produce for instance a A. This A 
is normally strongly suppressed because of the high mass, but now it delivers the baryon 
number N = 1 and is therefore enhanced. In order to balance the K°, the system has to 
produce an A and antibaryons are strongly suppressed. This causes the difference in the 
K° and K° canonical factors. 
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6.3 Heavy particles in the canonical ensemble 

All particles containing charm or bottom quarks are heavy and therefore suppressed by the 
Boltzmann factor. Thus the thermal production of heavy particles is very small^. In the 
canonical ensemble the canonical factor affects the particle production and depends on the 
mass of the lightest particle which has to be produced to balance the quantum numbers. 
In the case of charm particles, the system has to produce another charm particle and the 
lightest possibility is the with a mass of 1864 MeV. Therefore the canonical suppression 
is expected to be very high. This is shown in Figure [631 
In a neutral system the canonical factor for a is 

ZC(X) ZC(0) \ I I 

The exact value depends on the parameters (the given range is close to the parameters 
in e+e^ collisions), but in any case there is a strong suppression for charm particles in 
a neutral system by the canonical factor, in addition to the suppression of the Boltzmann 
factor caused by the high mass. 

However, as already discussed in the last section, the canonical factor can also enhance a 
particle yield if this is necessary for quantum number conservation. In a system with initial 



^The particle partition functions of the D'-* and are 0(10 ^) and 0(10 ^^), respectively. 
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charmness (C=l), like a c-jet in an e+e collision, the situation is exactly inverse to what 
is shown above 

£!™.£!^. 0,100,-0(1000) (6.10) 

In a charmed system all charm particles are strongly enhanced because they are needed to 
fulfill the initial charmness. This is also shown in Figure [631 and obviously the canonical 
factor plays a dominant role in the production of heavy particle^. 



6.4 Other quantities 

The grand-canonical partition function i^ 

InZ^^ = J2 I d'p e-"^'^-^^^ (6.11) 

where the Boltzmann approximation is applied. By using eqs. (12.221) . (|2.26l) and (|2.27l) it 
is possible to calculate the average energy, entropy and pressure of the system 



^GC ^ _ 



dp 



^t=0 J J 



S^^ = ^^^^^ = lnZ^^ + ^Y: (6-13) 

j 

,ac a(riiiZGC) g, 



The canonical partition function is (eq. (14.411) for simplicity in pure Maxwell-Boltzmann 
statistics) 



with the particle partition function 



'^In the case of bottom particles the canonical factor in a system with bottomness (-6=1) is w 0(10^"), 
because of the heavy mass of the bottom particles. 
^The general case /i 7^ is treated here. 
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In the canonical case, the energy, entropy and pressure are 



E 



c 



df3 



ZC{X) 



J l/ij=0 



S 



/ij=0 



dT 



dV 



E 



Z {X Xj) ^Qfj, 



T 



(6.15) 



(6.16) 



(6.17) 



At each quantity the canonical factor supplements the grand-canonical quantity, very sim- 
ilar to the situation for the multiplicities where 



inf) 



dlnZ 



c 



Xi=l 



(6.18) 



was found (see chapterS]). 



Chapter 7 

Measurement properties and the fit 
procedure 



In this chapter the fit procedure to describe the experimental particle spectra is presented. 
The contribution of strong and weak decays is discussed as well as the strangeness under- 
saturation parameter and the influence of the particle width. 



7.1 Particle width, strangeness under-saturation and the 
treatment of strong decays 

Within the calculation of the particle partition function, the particle width has to be taken 
into account. This is important, because the number of lightparticles coming from the 
decay of resonances is increased by the finite resonance widtlu. The finite widths of res- 
onances are taken into account by an additional integration over the particle mass in the 
multiplicity calculation, with a Breit-Wigner distribution as a weight 11711 



^1 = f9\\3 AT / ^"^7 \2^r2/A / ^ ^ ~ ''^^ 

{2Tchy Nbw J {m-mjY + TyAJ 

Mo 

where rrij is the nominal mass and Tj is the width of particle j. The energy is then calcu- 
lated for every value of m in the integration step. Here, Nbw is the normalization of the 
Breit-Wigner distribution such that the integral over the Breit-Wigner factor gives 1 . Mq 
is the threshold defined by the masses of the particles in the dominant decay channels. For 
instance, at the A++, there should be strength (Breit-Wigner amplitude) below a thresh- 
old defined by the sum of the masses of the dominant decay channels (proton and pion). 
Following the approach of ref.[6J, an additional parameter 7^ is introduced into the particle 



'The effect is decreasing with increasing temperature. 
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7 Measurement properties and the fit procedure 



partition function to account for a possible deviation of strange particle yields from their 
chemical equilibrium values (see section I3.3.3h. If a hadron contains Us strange valence 
quarks, its production is reduced by a factor 7"'! This parameter is also applied to neutral 
mesons such as t], t]', 0, cu, /2(1270) and /2(1525) according to the fraction of ss content in 
the meson itself. The relevant fraction is determined by using mixing formulas quoted in 
10. 

The multiplicity calculation basically proceeds in two steps. First, a primary hadron yield, 
(n^'*) is calculated using eqs. (14.241) and (|4.47l) . As a second step all resonances in the 
hadron gas which are unstable against strong decays are allowed to decay into lighter sta- 
ble hadrons, using appropriate branching ratios (Br) and multiplicities (M) for the decay 
j ^ h published by the PDG [5J. The abundances in the final state are thus determined by 



) = «') + $^(r2r')Br(j - h)M{j ^ h) (7.2) 



where the sum runs over all particles contained in the hadron gas and h refers to a hadron. 
For many particles the decays give significant contributions to the final yields. For example 
the measured 7r° yield consists of about 25% initially produced 7r°, whilst the rest is from 
(mostly strong) decays (for e~^e~ collisions at 91 GeV). 

A stable particle is a particle which does not give any decay contributions because it lives 
long enough to be detected itself, or at least forms a secondary vertex which can be sep- 
arated from the primary vertex, such that the original particle can be reconstructed. The 
proton, neutron and the three pion states are treated as stable particles. The weak decaying 
strange particles are discussed in the next section. 



7.2 Weak decays 



Particles which decay by the strong or electromagnetic interaction preserve the strangeness 
quantum number. The decay process for e.g. the A particle must violate this rule, since 
there is no lighter baryon which contains a strange quark, so the strange quark must be 
transformed into another quark in the decay process. That can only happen via the weak 
interaction where strangeness is not conserved, and this leads to a much longer lifetime. 



uds 
A 

uds 
A 



uud 
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udd 
n 
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+ 
+ 



du 

uu+dd 



> dominant A decays 



Strangeness: —1 7^ + 







^Please take into account that for the quantum statistical calculation with terms fc > 1, 7^"' has to be 
used. 
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Table 7.1: Decay rates of the weak decaying particles (within cr=10 cm). The mo- 
mentum spectra were extracted from the given references. The values 
with * are estimates, because no momentum spectrum was available. The 
lifetimes are taken from the PDG tables flU. 



The long lifetime of the A rises the question whether there are As or its decay products that 
are seen in the detector. The PDG data used in this analysis include decay products from 
all resonances with decay lengths cr < 10 cm [5]. Therefore approximately 30% of the As 
will survive and give no decay contributions. There is a similar behavior in the case of , 
Kg, , S^, S+ and n~ (and the corresponding antihyperons). This fact is accounted 

for in the present calculations. The exact values implemented in the model are shown in 
Table 17711 The importance of the effect can be seen at the K^, which is an almost stable 
particle (in the detector sense of stable), because of the long lifetime, whereas almost 100% 
of the S+ decay. 

7.3 The effect of uncertainties in the particle table 

The particle table is surely one of the biggest uncertainties of the model, because 

• Above 2 GeV the particle spectra is not fully measured. It is assumed that it increases 
exponentially 1(76117711781. 

• The decay channels, especially for heavy particles, are not well known. 

Fortunately the importance of the particles (for our results), decreases with increasing mass 
because the Boltzmann factor suppresses all heavy particles. Therefore it is acceptable to 
use an upper cut of the spectra with the assumption that all particles above this cut can be 
neglected. Nevertheless, in the correlated scheme, large enhancements for heavy particles 
can occur. Such results are connected to large uncertainties. 

To investigate the effect of an upper mass cut, the calculation was performed with different 
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Table 7.2: Particle yields at ^/s=9l GeV for different mass cuts. Please take into 
account that above 2 GeV the mass spectra is incomplete and the higher 
mass cut includes only the known states of [5j. 



mass cut values. The results are shown in Table 17.21 The effect seems to be small, but it is 
of the order of the experimental errors. 

The set of hadron species used in this work was updated compared to IItTI with the most 
recent information available from the PDG [|5l and includes all resonances listed there. 
Concerning the decays, it was ensured that the decay tables are symmetric for particles 
and antiparticles and the decay widths always add up to the total width, even if particu- 
lar channels are not measured. As the branching ratios are not well known for some of 
the high-mass states, the known Br of the nearest state with the same quantum numbers 
are used. The particle spectra is shown in Figure 17.11 and it is clearly visible that the ex- 
pected exponential increase does not hold for masses above 2 GeV. This is caused by the 
incompleteness of the measurements. 



7.4 The fit function 



After calculating the particle yields using the equations of chapter |4] and taking into ac- 
count all the influences discussed in this chapter, the particle yields can be compared to the 
measurements. All fit procedures are carried out by minimizing (least x^"fit) 



X 



(7.3) 
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Figure 7.1: Mass spectra of the used mesons (red) and baryons (blue). 



as a function of the three parameters T, V and 7^, taking account of the experimental 
uncertainties aj. aj contains only the errors of |l5l. There is no method applied to include 
the uncertainties in mass, width and branching ratios as additional systematic errors like 
done inilleolllll. 

The value {nj°^) depends on the used scheme (correlated or uncorrelated) as explained in 
section l44l The schemes used in this analysis are discussed in the next chapter. 



Chapter 8 

Results and discussions 



In this chapter the results of the data analysis of e^e^ collisions at ^/s=lO GeV, 29-35 GeV, 
91 GeV and 130-200 GeV are discussed. The model is based on the canonical partition 
function (in quantum statistics) with the conservation of five quantum numbers as derived 
in section l4".2.1[ The fits include all light particles measured as given in appendix|Dl Heavy 
particles will be treated separately in section [831 

8.1 The model 

The analysis is focusing on thermal production and is neglecting the issue of hard colli- 
sions. Therefore an uncorrected jet scheme is used, i.e. each jet is treated as a fireball 
with vanishing quantum numbers as fixed by the entrance channel. It is clear at this point 
that hadrons from jets with heavy quarks (c and b) will be greatly underestimated by the 
model because of the large Boltzmann suppression factors. In this approach the issue of 
equilibration is effectively addressed only for hadrons with light quarks (u, d, s). 
It is important to recognize that the measured yields of these hadrons contain the con- 
tribution from the e+e annihilation events into cc and bb. Heavy-quark production is 
indeed significant and is very precisely measured, in particular at the mass (^/s=9l.2 
GeV), where the measurements are well described by the Standard Model [5]. Hence, 
heavy-quark production is manifestly non-thermal in origin. Therefore two scenarios were 
considered: 

i. Fit to the data as measured. 

ii. Subtraction of the contribution originating from charm and bottom decays from the 
yields of hadrons carrying light quarks based on available data for charm and bottom 
hadron production (and their branching ratios) at 91 GeV. 

This approach differs from that used by Becattini et al. (see [80J for the most recent re- 
sults). There, e+e^ cc and e+e~ — > bb events are treated in a correlated jet scheme 
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Table 8.1: Comparison of particle yields obtained with the code used in this analysis 
(in quantum statistics (QS) and Boltzmann statistics (BS)) and with the 
THERMUS code [9 8 J for both initial production and for the final values 
(after strong and electromagnetic decays). The parameters listed corre- 
spond to the best fit at 91 GeV (r=157 MeV, V=32 fm^ and 7,=0.80). 
No contribution from weak decays is included. 



(a 2-jet initial state which carries quantum numbers C=±l or i?=±l). The fractions of 
hadronic events with 6 or c quarks, as well as those with lighter quarks, are thus external 
input values, unrelated with the Thermal Model (see Table II in ref. [80]). For the light 
quark sector, the difference between calculations with uncorrected and correlated quark- 
antiquark schemes are investigated later in this chapter. 

The second scenario, namely subtracting the contribution originating from charm and bot- 
tom decays, is investigated only at v^=91 GeV as a case study, as the measurements 
needed to allow the subtraction are complete only at this energy. 

8.2 Comparison to other calculations 

Before engaging in the present data analysis a comparison of results from the model de- 
scribed above with those obtained using the THERMUS code ll98l was performecfl. The 

'The THERMUS code is available publicly and can be downloaded from the web: 
|http : / /hep . phy . uct~ ac . za/THERMUS/SourceDownload . html| 
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8 Results and discussions 



results for a particular set of parameters are presented in Table 18.11 The deviation caused 
by the Boltzmann approximation is also shown in this Table. Here and in the following the 
quoted yields include the corresponding antiparticle states (i.e. the yield labeled 7r+ is ac- 
tually the sum of the yields of vr"*" and vr^) and are the sum over the two jets, as summarized 
and published by the PDG f5\. The agreement between the particle yields obtained with 
both codes is very good, lending strong support also to the numerical implementation of 
the methods used in this analysis. Note that, in the calculations used here, quantum statis- 
tics is employed, whereas in the THERMUS code the Boltzmann approximation is used 
for the canonical ensemble. This causes an error of about 3% for the final pion yielcH i.e. 
larger than the uncertainty in the data. As expected, the effect on all other hadron yields 
is smaller. If the code used in this analysis is restricted to Boltzmann approximation the 
yields, prior to strong decays, for mesons and non-strange baryons are in agreement with 
THERMUS results at the percent level, while for strange baryons the yields are system- 
atically higher by about 2% compared to those obtained with THERMUS. Inspecting the 
yields after strong decays, one notices a discrepancy between results from both calcula- 
tions and for most of the hadrons. The reason is a more complete set of hadron yields used 
in the calculations. 

The measured hadron yields summarized in Q include decay products from all resonances 
with decay lengths cr < 10 cm. The effect to weak decaying particles is appropriately taken 
into account in the present calculations as shown in Table 17.1 1 Because of weak decays the 
final yields are increased by 10% for charged pions (8% for 7r°) and 23% for protons. The 
magnitude of this contribution underlines its importance for the analysis of data. 
In the resonance gas model the results are determined by the basic thermal parameters, 
temperature T and volume V (the volume corresponding to one jet). 

8.2.1 Comparison to former results 

The comparison of the present results with those published in |l6l |9TJ [601 is rather unsatis- 
factory, with discrepancies exceeding 50% for many particle yields (see Table l8.2l) . For the 
parameters of the best fit of ref. [6] in the uncorrelated jet scheme (T=166.7 MeV, 1^=17.5 
fm'^, 7^=0.698) for ^/s=9l GeV, the 7r+ yield from the code used in this work is 6.87, 
while that published in [6J is 8.72, whereas for the ^7 baryon the corresponding yields are 
2.60-10"'' and 7.82-10"^, respectively. This is also The THERMUS results are included 
in Table 18.21 to show the good agreement with the results of this work. Each deviation 
between the results obtained in this work and the THERMUS results can be explained by 
a more complete set of particles. 
The model differences of this work to [6J are: 

1. Conservation of five quantum numbers neglects the conservation of electric 
charge). 



^The differences discussed in section |52] are from the initial yields. 
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T=166.7 MeV, ^=17.5 fm^, 7,=0.698 



particle 


results of [6] 


This work 


IHrLKMUa 




8.72 


6.33(6.87) 


5.55 


n 


9.83 


7.55(8.05) 


6.41 


K+ 


1.06 


0.653(0.653) 


0.607 


P° 


1.16 


0.948(0.948) 


0.878 


P 


0.484 


0.301(0.377) 


0.290 


0(1020) 


0.167 


0.0916(0.0916) 


0.0873 


A 


0.152 


0.120(0.130) 


0.0864 




0.0110 


0.00557(0.00560) 


0.00524 


Q 


0.000782 


0.000260(0.000260) 


0.000247 



Table 8.2: Model comparison, to THERMUS [98] and [6]. The values in brackets 
include weak decays, whereas THERMUS yields do not include weak 
decays. 



2. Quantum statistics is applied, whereas ^ uses Boltzmann statistics. 

3. The particle table includes all particles published by the PDG (2008) with a mass cut 
at almost 3 GeV, whereas in [i6J a mass cut at 1.7 GeV is used based on the data of 
1996. 

Unfortunately the differences of the values obtained in this work to the results of |l6l cannot 
be explained by the three points mentioned above and are not understood so far. 

8.3 Fit to the experimental data 

The resulting best fit to the data at the energy ^/s=9l GeV is shown in Fig. 18. 1[ The two 
cases of the fit, without and with the subtraction of the contribution from heavy quarks (the 
magnitude of this contribution is listed in Table [84]) . are shown. The overall behavior of 
the data follows an approximately exponential decrease of particle yield with increasing 
particle mass. Such a behavior is expected in the hadron resonance gas model due to the 
Boltzmann factors, thus indicating the presence of statistical features of hadron production 
in elementary collisions. 

The quantitative description of the data with the Statistical Model is, however, rather poor 
and certainly no improvement is visible for the case of subtracting charm and bottom con- 
tributions (scenario ii). The poor fit quality which is already visible in Fig. 18. II becomes 
striking when investigating Fig. 18. 2[ where for the four energies the difference A (in units 
of the experimental error) between the experimental data and the Statistical Model calcu- 
lations for the best fit values are shown. 
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Figure 8.1: Comparison between the best fit Thermal Model calculations and ex- 
perimental hadron multiplicities for e+e^ collisions at v^=91 GeV. The 
upper panel shows the case including in the data the contribution from 
heavy quarks, the lower panel is after subtraction of this contribution. 



A summary of the fit parameters obtained for the various data sets is presented in Table l8.3l 
The errors of the thermal parameters listed in Table l8.3l are the result of a minimization per- 
formed with MINUIT ||99l . interfaced with the used code. For the true size of the errors, 
see the discussion in connection to Fig. l8.3l later in this chapter. 

Typical values per degree of freedom lie between 5 and 20 and discrepancies between 
single data points and fit values larger than 5 standard deviations are not rare. Furthermore, 
there is no clear pattern observed: the fits are comparably poor for baryons and mesons, as 
well as for non-strange and strange hadrons. In particular, for all energies the yields of 
mesons and A and S hyperons are poorly reproduced. Large deviations are seen also for 
kaons. 

The case of the fit after subtraction of the decay contribution of charmed and bottom 
hadrons (scenario ii), explored at v^=91 GeV, is characterized by a larger value com- 
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Figure 8.2: Difference (in units of experimental error) between experimental data 
and thermal model fits at four energies. At ^/s=9l GeV both scenarios 
are shown, the uncorrected scheme (black point) and the results of the 
fit to data after subtraction of the heavy quark decay contribution (open 
circles). 



v/g[GeV] r[MeV] V [fm^] 7, x^/dof 

10 159±1.7 14±1.5 0.80±0.02 325/21 

29-35 160±1.7 18±1.4 0.96±0.03 105/18 

91 (all) 157±0.50 32±1 0.80±0.007 659/30 

91(-c,b) 169±0.50 18±1 0.68±0.01 781/30 

130-200 155±2.8 40±4.3 0.78±0.03 12.5/2 



Table 8.3: Values of fit parameters and values per degree of freedom for different 
energies. 
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68.0 


A 


0.343 


0.363 


0.383 


1.5 


S+(1385) 


0.0392 


0.0416 


0.0438 


~0 




0.0152 


0.0172 


0.0177 


~0 


Q 


0.000751 


0.000936 


0.000960 


~0 


X^/dof 


781/30 


762/30 


716/30 





Table 8.4: Calculated yields (corresponding to the fit of the subtracted data at 91 
GeV, T=169 MeV, 1^=18 fm^ and 7^=0.68) assuming vanishing or frac- 
tional baryon and charge quantum numbers. Also x^/dof of the fits for 
the three cases is shown. In the rightmost column the percentage in- 
crease of the yields due to the contribution of charm and bottom events 
is shown, as calculated from the experimental data. 



pared to the overall fit (scenario i). The extracted fit parameters differ in the two cases. In 
particular, T is higher and V is smaller for scenario ii) compared to scenario i), following 
the (T, V) correlation shown in Fig. 18.31 

It should be mentioned that while the general agreement among the four LEP experiments 
is excellent, the measured yields of E* hyperons differ by more than 70%. Excluding the 
S* hyperons would cause an increase of 7s at ^/s=9l GeV, because the S* yields calcu- 
lated in the model are overestimated (see Fig. 18.11) . This would also slightly improve the 
situation for the Q and A multiplicities which are much higher than those predicted by the 
model. However, excluding those particles from the fit does not result in a significant im- 
provement of the values, as is discussed below. 

Now the difference between calculations with uncorrelated and correlated quark-antiquark 
(u,d,s) schemes will be discussed. In this context one should note that in the correlated 
jet scheme there is an open issue concerning the treatment of electric charge and baryon 
numbers of the quarks. A calculation for vanishing charge and baryon number and for 
fractional quark quantum numbers is shown in Table l8.4l The results in these 2 cases differ 
little for all hadrons, with the exception of protons, for which a difference of 11% is ob- 
served. However, the quality of the fit to the data (charm- subtracted data at 91 GeV) is not 
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Figure 8.3: contour lines in temperature and volume space for the overall fit (left 
panel) and after subtraction of charm and bottom decays (right panel) 
for ^=91 GeV. The contours correspond to xLn+10%, +20%, +50%, 
+100%. The best fit values are indicated by the blue crosses. 



significantly improved. 

A further difficulty is visible if one inspects contour lines as shown in Fig. l8.3l in (T,V) 
space for both fit scenarios at 91 GeV. One notices in this Figure a strong anti-correlation 
between the fit parameters which is also present in the (T,7s) space (shown in appendix 
IC.2I) . Closer inspection reveals, in addition, a series of local minima which indicates the 
difficulty in the determination of the fit parameters. Such local minima are typical for poor 
fits and imply that the true uncertainties in the fit parameters are likely much larger than 
the values obtained from the standard fit procedure [|99l employed here. 
Despite these caveats about fit quality and uncertainties it is noteworthy that the tempera- 
ture parameters obtained from most data sets are close to 160 MeV and nearly independent 
of energy, similar to results of previous investigations. In contrast, the volume increases 
with the center of mass energy, reaching a value at 91 GeV which is larger by almost a fac- 
tor of 2, compared to results reported in O, but somewhat smaller than the value obtained 
in [|9T1 1601 . These differences may also be connected with the local minima observed in 
the surface discussed above. The values obtained for the strangeness under- saturation 
parameter 7^ range between 0.96 and 0.68 and exhibit no clear trend with energy. 
In the following the 91 GeV fit will be investigated in detail, as the data set at that energy 
contains the largest number of measured hadron yields. Since there is no statistically sig- 
nificant difference between the two fit scenarios, scenario i) is considered, that is, the fit to 
data without subtraction the charm and bottom decay contributions. To check whether the 
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high x^/dof values are caused by discrepancies for a few particular particles, the four worst 
hadron species (the three S*'s and 0) are excluded. Still a high x^/dof=206/26 (for T=159 
MeV, V=2% fm^ 7^=0.86) remains. 

The data from LEP comprise many particle species and their yields typically are measured 
with accuracies of a few percent. Clearly the precision of the experimental data set provides 
a stringent test of the Statistical Model. To provide a quantitative estimate at which accu- 
racy level the Statistical Model breaks down and to allow comparison with the situation 
encountered for RHIC data a fit at y/s=9l GeV is performed to those hadron yields which 
were used for the fit of central nucleus-nucleus collision data at y^J/v7v=200 GeV [fTTI . 
namely 7r+, n^, K+, K*+, K*°, p. A, E', fi", and 0. 

The resulting thermal parameters are: T=162±2 MeV V=26±2 fm^, and 7s=0.78±0.02, 
with a still poor x^/dof of 263/8. Excluding the 0-meson from the fit yields T=160±2 MeV, 
V=26±2 fm^ and 7s=0.90±0.02, with x^/dof=57/7. A more direct comparison between 
thermal fits for heavy ion and e~^e~ data can be obtained by assigning the uncertainties 
of the RHIC data to the corresponding LEP yields (namely uncertainties of the order of 
13%). Constraining, as for the RHIC data, the fit parameters to T and V, i.e. setting 7s=l 
yields then r=157±2 MeV and V"=26±2 fm^, with x^/dof=38/9. A reasonable fit can 
only be obtained by also letting 7^ vary freely, with resulting parameters T=166±2 MeV, 
\^=22±2fm^ and 7s=0.76±0.02, with x^/dof=l 1.8/8. These exercises demonstrate that a 
Statistical Model description of e+e" data fails badly without the introduction of the non- 
equilibrium parameter 7^. Even using 7s the statistical description breaks down completely 
at an accuracy level for the data of below 10%. 



8.4 Energy content 

Another noteworthy difference between fireballs in e+e^ and nucleus-nucleus collisions 
is their energy content. For its determination the energy density e in the hadronic gas is 
computed to yield the thermal energy content E=eV of one jet at chemical decoupling. 
The energy can be derived by applying eq. (|2.22l) to eq. (14.411) 

d\nZ 

^ Z{X) ^ Z{X) ' 

with 

where j runs over all particles and h runs over the light bosons. The equation shows, that 
the canonical factor will decrease the energy as in the case of the particle yields. 
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Quantum 


Boltzmann 


Grand 


[GeV] 


statistics [MeV] 


statistics [MeV] 


canonical [MeV] 


10 


3183.07 (60.1%) 


3169.46 (59.8%) 


5298.64 


29-35 


5380.70 (69.9%) 


5354.05 (69.6%) 


7697.03 


91 


8773.02 (80.1%) 


8713.98 (79.6%) 


10951.16 


130-200 


10159.63 (82.9%) 


10085.10(82.3%) 


12256.37 



Table 8.5: Energy content of one fireball for the four different energies, calculated in 
quantum statistics, Boltzmann statistics and grand-canonical. The com- 
parison to the grand-canonical result are given in %. 



The results for the parameters reported in Table 18.31 are shown in Table 18. 5[ The size 
of the compression because of the canonical factor decreases with increasing energy, like 
one would expect. It is also possible to see the small but visible effect caused by the un- 
derestimation of Boltzmann approximation in comparison to the correct quantum statistic 
calculation. Note that the thermal energy within each jet is only a small fraction of -\/i/2, 
(e.g. about 19% at 91 GeV). Apparently, in the thermal interpretation of e+e~ collisions, 
most of the cm. energy is not available for particle production. This is in strong contrast 
with results for nucleus-nucleus collisions. For this purpose central collision events for 20 
and 40 GeV/nucleon Pb-Pb collisions [71] are analyzed, where it makes sense to consider 
data integrated over the full phase space. In these cases the energy content in the fireball 
amounts to 61% and 63% of the total cm. energy at 20 and 40 AGeV, respectively, imply- 
ing that most of the total cm. energy in a nucleus-nucleus collision is thermal, with the 
remaining non-thermal fraction likely to be due to collective flow. 

These differences are not consistent with the finding that particle production in e+e~, 
pp and nucleus-nucleus collisions is universally governed by the available cm. energy 
[[Tool [Ton [T021 . 

8.5 Heavy flavored particles 

The thermal production of particles containing charm or bottom quarks is very small, but 
as already shown in section [63] there is a strong enhancement of charm particles in c-jets 
and bottom particles in b-jets. The question is, whether a correlated jet scheme, with the 
appropriate treatment of quantum numbers in the corresponding jets is able to describe the 
measured heavy particle yields. 

Figure 18.41 and 18.51 show the measured charm and bottom yields together with the calcu- 
lated values in a correlated (thermal jet fragmentation) and uncorrected (thermal model) 
jet scheme. The uncorrelated scheme corresponds to the pure thermal production. The 
parameters are T=169 MeV, V=IS fm^ and 7s=0.68 as found for the c, b subtracted fit to 
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Figure 8.4: Model predicted values for charged particles, compared to the LEP mea- 
surements at 91 GeV. The black lines corresponds to the uncorrelated jet 
scheme and the red lines (thermal jet fragmentation) corresponds to the 
correlated jet scheme. 



the light particles. In contrary to the uncorrelated scheme (black lines) the correlated jet 
scheme (red lines) is able to describe the measured data with the exception of the T yield 
in Figure 18. 5[ which can only be produced thermally even in the correlated jet scheme, 
because it has no net bottomness. The J/\l/ yield can be described because the main con- 
tribution of it comes from bottom decays. The x^/dof of the four bottom particles (without 
the T) is 0.2/4 and of the charm particles it is 76/10. Nevertheless it seems, that it is almost 
independent of the parameters what can be seen in Figure [831 

It was already discussed, that there is the necessity of an additional parameter 7^ for par- 
ticles carrying strangeness caused by an incomplete equilibrium. The reason is the high 
mass of the strange quark and therewith the Boltzmann suppression. The mass of the 
charm quark is much higher than the strange quark mass, thus there is no expectation of 
an equilibrium of charm particles. So why is it possible to describe the charm and bottom 
particle yields without any under- saturation factor as used for strange particles? 
The thermal production of heavy quark is very low (black lines in Figure [84l and 18. 51) . The 
charmness and bottomness production takes place in hard collisions. This is taken into ac- 
count by the parameters Rq. The sum over all initial charm particles gives exactly Rc=0. 17 
(or Rb=0.22 for bottom particles). This can be seen by investigating the charm particle 
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Figure 8.5: Model predicted values for bottom particles, compared to the LEP mea- 
surements at 91 GeV. The thermal model corresponds to the uncorre- 
lated jet scheme, thermal jet fragmentation is the results of the corre- 
lated jet scheme. 



yields using equation (14.381) and (14.131) which is in the case of a c-jet 

l^^charm\^ j N—Nj,S—Sj,Q—Qj 

(and similar for bottom particles). 

This is not true for u, d and s quark, because the thermal production of this particles is much 
bigger than the production in hard collisions (In the case of strange particles there is indeed 
a conflict between the parameter 7^ and Rs). The production of strange and charm particles 
in the different jets is shown in Table 18. 6[ It can be seen, that even in the antistrange jet 
there is a production of strange quarks, because strangeness is produced thermally. One 
could ask why is there a difference of strangeness production in e.g. the u and u jets. This 
differences are caused by the electric charge. The values in Table 18.61 are the sums over all 
initial produced strange (S'=-l) particle yields. This is e.g. in case of kaons the and the 
K°. and it is clear that there are more in an w-jet than in an w-jet. Similar for the other 
jets, whereas the difference in the s-jets is obviously caused by the strangeness quantum 
number. It is again visible, that the strangeness quantum number plays a much bigger role 
than the electric charge quantum number, what was already discussed in chapter [6] and 
shown in Figure IGTI 
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jet 


strangeness (S=-l) 


charmness (C 


M-jet: 


0.126 


2.184-10-6 


M-jet: 


0.127 


2.184-10-6 


djet: 


0.163 


2.809-10-6 


rf-jet: 


0.162 


2.809-10-6 


s-jet: 


0.322 


3.06-10-6 


s-jet: 


0.108 


3.06-10-6 


c-jet: 


0.144 


0.17 


c-jet: 


0.145 


6.776-10-=^° 


6-jet: 


0.191 


3.036-10-6 


6-jet: 


0.189 


3.036-10-6 



Table 8.6: Sum over all initial strange/charm particle yields in the different jets. The 
values are already multiplied wit Rg (see eq. (14.591) ). 
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In contrary to strange particles thermal production of charm particles is very small and the 
final charm yields are dominated by the production in c-jets. The overall charmness in such 
a jet is defined by the parameter Rc. 



Chapter 9 

Summary and Conclusion 



Within this thesis a numerically calculable quantum statistical canonical partition function 
taking into account the conservation of four and five quantum numbers was derived |I3. 
Based on these functions a computer program was implemented to calculate the hadron 
production in particle collisions [j4]|. Furthermore the canonical factor with its special fea- 
tures was investigated and explained in detail. The comprehensive new set of measured 
yields of hadrons with light quarks {u, d, s) in e+e" collisions was analyzed in a range 
of collision energies from 10 GeV up to 200 GeV Q. At ^/s=9l GeV, two scenarios 
were considered, namely fitting data without and with subtraction of the decay products 
of charmed and bottom hadrons. The two thermodynamical parameters temperature T and 
volume V, and a strangeness under- saturation factor 7^ were obtained from a minimiza- 
tion procedure. While, as in previous investigations ll6l|9T]|60l, the resulting temperature 
value was close to 160 MeV, independent of energy, the overall description of the high- 
precision LEP data is rather poor, independent of whether heavy quark contributions are 
subtracted or not. The x^/dof values are larger than 5 for all fits and call into question the 
validity of the thermodynamical approach for these data. This conclusion still holds even 
if the analysis is restricted to the same set of hadrons which were analyzed in the context 
of thermal model fits to Au-Au collision data from the RHIC accelerator. 
The apparent statistical fingerprint visible in the LEP data and first observed in (6] breaks 
down at an accuracy level of about 10%. Even at that level the e+e^ data cannot at all be 
described without the explicit assumption of strangeness under-saturation, implying that 
hadron production in e+e~ originates from a state which is quite far from true thermo- 
dynamic equilibrium. This conclusion is further supported by the observation that the 
corresponding fireball volume contains only a small fraction of the overall center of mass 
energy, in strong contrast to the situation in nucleus-nucleus collisions, where the fireball 
volume contains the majority fraction of the available energy. 

Furthermore it was observed, that the implementation of fireballs with charm and bottom 
quantum numbers can be used to describe the heavy particle production in e+e~ collisions 
[|2| determined by the new input values Rg (except of the T yield). This agrees with the 
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observation of [l6l[80l|. However, the contours plot showed that a wide parameter range of 
volume and temperature is able to explain the measured data. Nevertheless it seems that 
the production of this particles is indeed controlled by the particle phase space and the 
canonical factor. 



Appendix A 

Mathematical and numerical 
approximations 



The following integral has to be solved (I4.42|) 



k{2nhy 



because the integral depends only on the absolute value of p, it simplifies by changing into 
spherical coordinates 



k 



dr r^e t ^ (A.2) 



This integral can be solved numerically. 

It is also possible to use the following relation 



with the modified Bessel-function of order 2. This relation obeys the following limiting 

cases 

m\2 /m\ f2 for m^O 

and equation (|A.2I) becomes 



K2[-) ^ \ ^ (A.4) 

TJ \T J /^g-T for m>T 
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A.l Calculation comparison 



particle Bessel-way Integral-way 
(section |42lB (section |4]3]) 



7r± 


3.186 


3.185 




3.629 


3.634 


V 


0.510 


0.510 




0.779 


0.779 


0(1020) 


0.1243 


0.1243 


P 


0.1288 


0.1285 


A 


0.03602 


0.03593 


E+ 


0.02230 


0.02225 




0.005667 


0.005690 


Q- 


0.0005432 


0.0005494 



Table A.l: Comparison of the two ways of calculation explained in chapter HI The 
shown yields are the initial ones, what means that there is no treatment of 
any decay. The deviations are caused by numerical approximations (uncor- 
related scheme). 



Appendix B 

Quantum statistical calculation (for 
bosons and fermions) 



The quantum statistical canonical partition function is (14.391) 



f)2 4_T-(-) 2 



as already shown in section 14.2.11 where the correct quantum statistical calculation was 
implemented for bosons only (14.401) 

oo ^ 
k=l 

For fermions the following expansion has to be used 

CO ^ 

ln(l + a;) = ^(-l)^+i^ (B.l) 

k=l 

By including the whole series for both, fermions and bosons, to the partition function one 
gets 



-ikxj<j) 



j d'^e'^^hxp 1^ f;(±)'=+^zje-^^^^^^j (B.2) 



(27r) 



j k=i 
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where the + is for bosons and the — for fermions. The particle partition functions are 

defined like in (14.421) . 

The particle multiplicities are 
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(B.3) 



To get a numerically calculable term for the integrals in (IB. II) one has to apply the same 
manipulations as discussed in section lT.2.11 The result will be 
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with 



ui = N + nl + n2 + n3 + n5 + n6 + n7 + n8 + n9 + nlO 
+ ^ mrim + ^ nrin + ^ ono + ^ gn^ 
+ ^ rn^ + ^ sn^ + ^ trit + ^ ura^ 

1^2 = "S* + n4 — — ^6 — — 2n8 — 2r29 — 3nl0 
+ ^hnh + Y^ Irii - ^pup - ^quq - 
— 2sns — 2t?2t — Swriu 



(B.5) 



(B.6) 



80 B Quantum statistical calculation (for bosons and fermions) 

1^3 = Q + nl — n2 + 2n3 + nA + nQ — n7 — n9 — nlO 
+ krik + Irii + mrim — nUn 

+ 2ono + qUq — rn^ — t^t — (B.7) 
Inserting this integral in equation (14.471) gives the final partition function. 
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Appendix C 
Fit results 

C.l Particle yields 



particle v/i=10GeV v/i=29-35 GeV v/i=91 GeV v/i= 1 30-200 GeV 



7r± 


6.91 


10.97 


18.08 


21.07 




4.32 


6.60 


10.47 


12.10 




0.627 


1.36 


2.04 


2.36 




0.593 


1.30 


1.96 


2.28 


V 


0.501 


0.796 


1.11 


1.25 


P^(770) 


0.739 


1.13 


1.89 


2.18 


p°(770) 


0.518 


0.726 


1.13 


1.29 


/s:*±(892) 


0.175 


0.402 


0.604 


0.694 


i\:*°(892) 


0.172 


0.395 


0.594 


0.682 


P 


0.192 


0.459 


0.895 


1.05 


0(1020) 


0.0646 


0.125 


0.133 


0.142 


A 


0.0685 


0.213 


0.354 


0.409 


S+ 


0.0174 


0.0529 


0.0901 


0.104 


S+(1385) 


0.00746 


0.0228 


0.0391 


0.0452 




0.00277 


0.0133 


0.0201 


0.0235 


S0(1530) 


0.00106 


0.00509 


0.00762 


0.00882 


n- 


0.000105 


0.000879 


0.00124 


0.00147 



Table C.l: A random set of calculated particle yields (particle + antiparticle in two 
jets) for the best fit values (see Table 18.31) at the four different energies 
used in the analysis. The uncorrelated jet scheme is employed. 
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C.2 Plots 
C.2.1 10 GeV 



C Fit results 
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Figure C.l: Comparison between thermal model predictions and experimental par- 
ticle multiplicities for e+e~ collisions at ^/s=lO GeV. 




Figure C.2: contour lines in temperature and volume (left panel) and temperature 
and 7s (right panel) space for ^/s=\0 GeV. The contours correspond to 
xLn left: +10%, +20%, +50%, +100% right: +1%, +2%, +5%, 
+10%, +20%. The best fit values are indicated by the blue crosses. 
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Figure C.3: contour lines in volume and 7^ space for v^=10 GeV. The contours 
correspond to xLn + 5%, +10%, +20%, +30%, +50%. The best fit 
value is indicated by the blue cross. 
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Figure C.4: Comparison between thermal model predictions and experimental par- 
ticle multiplicities for e+e^ collisions at y^=29-35 GeV. 
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C Fit results 




Figure C.6: contour lines in volume and 7^ space for yi=29-35 GeV. The con- 
tours correspond to xLn + 10%, +20%, +50%, +100%, +200%. The 
best fit value is indicated by the blue cross. 
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C.2.3 91 GeV (uncorrected) 
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Figure C.7: Comparison between thermal model predictions and experimental par- 
ticle multiplicities for e+e~ collisions at \/i=91 GeV. 
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C.8: contour lines in temperature and volume (left panel) and temperature 
and 7s (right panel) space for ^/s=9l GeV. The contours correspond to 
left xLn left: +10%, +20%, +50%, +100%, right: +1%, +2%, +5%, 
+10%, +20%. The best fit values are indicated by the blue crosses. 
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C Fit results 
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Figure C.9: contour lines in volume and 7^ space for yi=91 GeV. The contours 
correspond to xLn + 5%, +10%, +20%, +30%, +50%. The best fit 
value is indicated by the blue cross. 
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Figure C. 10: Comparison between thermal model predictions and experimental par- 
ticle multiplicities for e^e" collisions at -\/s=91 GeV. 
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V (fm^) "^s 
Figure C.ll: contour lines in temperature and volume (left panel) and temper- 



ature and 7s (right panel) space for ^/s=9l GeV. The contours corre- 
spond to xLn + 10%, +20%, +50%, +100%, +200%. The best fit 
values are indicated by the blue crosses. 




Figure C. 12: contour lines in volume and 7^ space for ^/s=9l GeV. The contours 
correspond to xLn + 5%, +10%, +20%, +30%, +50%. The best fit 
value is indicated by the blue cross. 
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C Fit results 



C.2.5 130-200 GeV 




Figure C.14: contour lines in temperature and volume (left panel) and temper- 
ature and 7s (right panel) space for yi= 130-200 GeV. The contours 
correspond to xLn + 10%, +20%, +50%, +100%, +200%. The best 
fit values are indicated by the blue crosses. 
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Figure C.15: contour lines in volume and 7^ space for y^=130-200 GeV. The 
contours correspond to xLn + 10%, +20%, +50%, +100%, +200%. 
The best fit value is indicated by the blue cross. 



Appendix D 



Data source [SQ 



D.l Heavy particles 



particle 


w 10 GeV 


Vi=29-35 GeV 


Vi=91 GeV 


Vs=130-200 GeV 


Pseudoscalar mesons: 








D+ 


0.194 ±0.019 


0.17 ±0.03 


0.175 ±0.016 




D° 


0.446 ± 0.032 


0.45 ± 0.07 


0.454 ± 0.030 




Dt 


0.063 ± 0.014 


0.45 ± 0.20 


0.131 ±0.021 




BW 






0.343 ± 0.032 




Bt 






0.057 ±0.013 




Vector mesons: 










D*(2010)+ 


0.177 ±0.022 


0.43 ± 0.07 


0.1937 ±0.0057 




D*(2007)'' 


0.168 ±0.019 


0.27 ±0.11 






i:'*(2112+) 


0.048 ± 0.014 




0.101 ±0.048 




5*(2) 






0.288 ± 0.026 




J/*(15) 0.0005 ± 0.00005 




0.0052 ± 0.0004 










0.0023 ± 0.0004 




t\is) 






0.00014 ± 0.00007 




Pseudovector mesons: 








Xci(3510) 






0.0041 ±0.0011 




Tensor mesons: 








5**(2) 






0.118 ±0.024 










0.0052 ±0.0011 










0.0083 ±0.0031 




Baryons: 










A+ 


0.074 ±0.031 


0.11 ±0.05 


0.078 ±0.017 




K 






0.031 ±0.016 






0.014 ± 0.007 









Table D.l: Average heavy particle multiplicities per hadronic e+e^ annihilation event at y/s « 10, 
29-35, 91, and 130-200 GeV. The rates given include decay products from resonances 
with cr < 10 cm, and include the corresponding anti-particle state 15J. (1) Contains 5+ 
and B°, (2) Any charged state (i.e. B*+, and B^) 
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D.2 Light particles 

D.2 Light particles 



91 



particle 


^/i w 10 GeV 


V^=29-35 GeV 


yG=91 GeV 


Vi=130-200 GeV 


Pseudoscalar mesons: 








^± 


6.6 ± 0.2 


10.3 ± 0.4 


17.02 ±0.19 


21.24 ±0.39 




3.2 ±0.3 


5.83 ±0.28 


9.42 ± 0.32 


— 




0.90 ± 0.04 


1.48 ±0.09 


2.228 ± 0.059 


2.82 ±0.19 


K° 


0.91 ±0.05 


1.48 ±0.07 


2.049 ± 0.026 


2.10 ±0.12 


V 


0.20 ± 0.04 


0.61 ±0.07 


1.049 ±0.080 


— 


77'(958) 


0.03 ± 0.01 


0.26 ±0.10 


0.152 ±0.020 


— 


Scalar mesons: 








/o(980) 


0.024 ± 0.006 


0.05 ± 0.02 


0.146 ±0.012 


— 






— 


0.27 ±0.11 


— 


Vector mesons: 








p(770)" 


0.35 ± 0.04 


0.81 ±0.08 


1.231 ±0.098 


— 


p(770)± 




— 


2.049 ± 0.026 


— 


w(782) 


0.30 ± 0.08 


— 


1.016 ±0.065 


— 


K*{892)^ 


0.27 ± 0.03 


0.64 ± 0.05 


0.715 ±0.059 


— 


K*{892)" 


0.29 ± 0.03 


0.56 ± 0.06 


0.738 ± 0.024 


— 


0(1020) 


0.044 ± 0.003 


0.085 ±0.011 


0.0963 ± 0.0032 


— 


Pseudovector mesons: 








/i(1285) 




— 


0.165 ±0.051 


— 


/i(1420) 




— 


0.056 ± 0.012 


— 


Tensor mesons: 








/2(1270) 


0.09 ± 0.02 


0.14 ±0.04 


0.166 ±0.020 


— 


/2(1525) 




— 


0.012 ± 0.006 


— 






0.09 ± 0.03 


— 


— 






0.12 ±0.06 


0.084 ± 0.022 


— 


Baryons: 










P 


0.253 ± 0.016 


0.640 ± 0.050 


1.050 ±0.032 


1.41 ±0.18 


A 


0.080 ± 0.007 


0.205 ± 0.010 


0.3915 ± 0.0065 


0.39 ± 0.03 


EO 


0.023 ± 0.008 


— 


0.076 ±0.011 


— 


E- 




— 


0.081 ±0.010 


— 


E+ 






0.107 ± 0.011 




E± 






0.174 ±0.009 






0.0059 ± 0.0007 


0.0176 ± 0.0027 


0.0258 ± 0.001 




A(1232)++ 


0.040 ±0.010 




0.085 ± 0.014 




S(1385)- 


0.006 ± 0.002 


0.017 ± 0.004 


0.0240 ± 0.0017 




S(1385)+ 


0.005 ± 0.001 


0.017 ± 0.004 


0.0239 ± 0.0015 




S(1385)± 


0.0106 ± 0.0020 


0.033 ± 0.008 


0.0462 ± 0.0028 




S(1530)0 


0.0015 ± 0.0006 




0.0055 ± 0.0005 




Q- 


0.0007 ± 0.0004 


0.014 ±0.007 


0.0016 ± 0.0003 




A(1530) 


0.008 ± 0.002 




0.0222 ± 0.0027 





Table D.2: Average hadron multiplicities per hadronic e+e^ annihilation event at y/s wlO, 29-35, 
91, and 130-200 GeV. The rates given include decay products from resonances with 
CT < 10 cm, and include the corresponding anti-particle state |5 1. 
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